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Abstract 

The effect of weak impurity disorder on flux lattices at equilibrium is stud- 
ied quantitatively in the absence of free dislocations using both the Gaus- 
sian variational method and the renormalization group. Our results for the 
mean square relative displacements B{x) = {u{x) — u{0))'^ clarify the nature 
of the crossovers with distance. We find three regimes: (i) a short distance 
regime ("Larkin regime") where elasticity holds (ii) an intermediate regime 
("Random Manifold") where vortices are pinned independently (iii) a large 
distance, quasi-ordered regime where the periodicity of the lattice becomes 
important. In the last regime we find universal logarithmic growth of dis- 
placements for 2 < d < 4: B{x) ~ j4dlog|x| and persistence of algebraic 
quasi-long range translational order. The functional renormalization group 
to 0(e = 4 — d) and the variational method, when they can be compared, 
agree within 10% on the value of A^. In d = 3 we compute the function de- 
scribing the crossover between the three regimes. We discuss the observable 
signature of this crossover in decoration experiments and in neutron diffrac- 
tion experiments on flux lattices. Qualitative arguments are given suggesting 
the existence for weak disorder in d = 3 of a " Bragg glass " phase without 
free dislocations and with algebraically divergent Bragg peaks. In d = 1 -|- 1 
both the variational method and the Cardy-Ostlund renormalization group 
predict a glassy state below the same transition temperature T = Tc, but with 
different B{x) behaviors. Applications to d = 2 + systems and experiments 
on magnetic bubbles are discussed. 
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I. INTRODUCTION 



The interest in the pinning of the Abrikosov vortex lattice by impurities was revived 
recently with the discovery of high-Tc superconductors. Impurity disorder conflicts with 
the long range translational order of the flux lattice and some glassy state is generally 
expected to appear. Understanding the precise nature of this new thermodynamic state 
and how it depends on the type of disorder existing in the system is very important for the 
determination of the transport properties of these materials, such as critical currents and 
I-V characteristicdi^. This problem, however, is only one aspect of the more fundamental 
and broader question of the effect of quenched impurities on any translationally ordered 
structure, such as a crystal. This question arises in a large number of physical systems under 
current active experimental study. Examples are charge density wavesi, Wigner crystalsi^ii 
magnetic bubble^i'0. Josephson junctionsElS, the surface of crystals with quenched bulk 
or substrate disordertZI, domain walls in incommensurate solidH. All these systems have in 
common a perfectly ordered underlying structure modified by elastic distortions and possibly 
by topological defects such as dislocations, due to temperature or disorder. The effect of 
thermal fluctuations alone on three-dimensional and especially on two-dimensional structures 
is by now well understood, and it was shown that topological defects are not important in 
the low temperature solid phased. Much less is known however on the additional effects of 
quenched disorder. In particular, the important question of precisely how quenched disorder 
destroys the translational long range order of the lattice is far from being elucidated. If 
disorder is strong, the underlying order is a priori destroyed at every scale and an analytical 
description of the problem starting from the Abrikosov lattice is difficult. One then has 
to use a more macroscopic approach based on phenomenological models such as the gauge 
glass mo dell'S'il. The success of these approaches then crucially depends on whether these 
effective models are indeed a good representation of the system at large scale, a largely 
uncontrolled assumption. If disorder is weak enough, however, one expects the perfect 
lattice to survive at short scales. Thus a natural first step for a theoretical description is 
to neglect dislocations and to treat the simpler problem of an elastic medium submitted 
to weak impurities. In that case one can consider a gaussian random potential created by 
many weak impurities with short-range correlations. In this paper we will focus on point-like 
disorder. The problem of correlated disorder, such as columnar or twin boundary pinning in 
superconductors!'!, which is also relevant for any type of quantum problems can be treated 



by similar methods, and is examined in details in Ref. 22 



This simpler problem of an elastic lattice in presence of weak disorder is already quite 
non trivial. Despite several attempts, its physics has not been completely understood. An 
important quantity, which measures how fast translational order decays, is the translational 
correlation function Cxoif) = {e^^o-{u{^)-'u.{o))''^ ^ where u{r) is the displacement from the 
perfect lattice and Kq one of the first reciprocal lattice vector. We denote by () and — , 
the thermodynamic average and the disorder average respectively. Ckq (^) can be extracted 
from the Fourier transform of the density-density correlation function at wavevectors near 
q = Kq, or directly measured by imaging the deformed lattice, and is thus a quantity easily 
accessible in experiments. 

A first calculation of Cit(,(r) was performed by Larkinll using a model in which weak ran- 
dom forces act independently on each vortex. These forces are correlated over a small length 
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^0, of the order of the superconducting coherence length. This model predicts that weak 
disorder destroys translational order below four dimensionsil and Cko{^) ~ exp(— r"*"*^). 
The destruction of the long range order in this simple gaussian model can be understood 
easily from the standard Imry-Ma argument. To accommodate the random forces a region 
of size R will undergo relative deformations of order u. The cost in elastic energy is 
while the gain in potential energy is uR'^/'^. The optimal is thus u = leading to 

the above decay of Cifp(r). Using similar arguments in the presence of an external Lorentz 
force, Larkin and OvchinikovS constructed a theory of collective pinning of the flux lattice. 
In this theory the critical current is determined from the length scale Rc at which relative 
displacements are or order u ^ C,o- This theory was very successful in describing conven- 
tional superconductors. However, a need to reconsider this theory was prompted by high Tc 
superconductors where the flux lattice is usually probed at larger scales. It turns out that 
the Larkin model, while it is useful for estimating critical currents, cannot be used to study 
large scale quantities such as translational order. 

In fact the purely gaussian model with random forces, and the resulting linear elasticity, 
becomes inadequate beyond the Larkin-Ovchinikov length Rc. It has only one trivial equi- 
librium state and responds linearly to external force. It is thus too simple to approximate 
correctly the full non linear problem and grossly overestimates the effect of disorder. At 
larger scales the lattice starts behaving collectively as an elastic manifold in a random po- 
tential with many metastable states, thus the exponential decay of Cxoir) in d = 3 cannot 
hold beyond r > Rc- Using known results on the so-called "random manifold" problem, 
Feigelman et al.i showed that the system presents glassy behavior and computed transport 
properties. This was also pointed out by Bouchaud, Mezard and Yedidia (BMY)c3'Eaj who 
used the Gaussian Variational Method (GMV) to study this problem, and found a power-law 
roughening of the lattice with stretched exponential decay of C/^p(r). 

However, the periodicity of the lattice was not properly taken into account in all the above 
works. The periodicity has important consequences for the behavior of correlation functions 
at large scales. Indeed, it was suggested with the use of qualitative Flory arguments, that 
periodicity leads to logarithmic roughening@, rather than a power law. 

In this paper we develop a quantitative description of the static properties of a lattice in 
the presence of disorder. A short account of some of the results of this paper were presented 
in a recent letterii. We take into account both the existence of many metastable states and 
the periodicity of the lattice. One of the difficulties in the physics of this problem is that the 
disorder varies at a much shorter length scale than the lattice spacing. As a consequence 
the elastic limit has to be taken with some care. Indeed in this limit the displacement varies 
slowly, but the density still consists in a series of peaks. To couple the density to the random 
potential it is thus important to distinguish between its slowly varying components Ag=o 
and its Fourier component close to the periodicity of the lattice. This separation of 
harmonics exposes clearly the physics and allows to treat all the regimes in length scales in a 
simple way. To study the resulting model, we mainly use the Gaussian Variational Method, 
developed to study manifold in random medial. We also use the Renormalization Group 
(RG) close to d = 4 dimensions and in d = 2 dimensions. Comparison of the two methods 
provides a confirmation of the accuracy of the GVM. 

One the main results of the present study, which is somewhat surprising in view of 
conventional wisdom based on Larkin's original calculation, is that quasi-long range order 
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survives in the system. This means that CKoir) decays as a power law at large distance. 
Such a property for disordered lattice is similar to the quasi-order found for clean two- 
dimensional solids. This state however has the peculiar property of being a glass with many 
metastable states, and at the same time show Bragg peaks as a solid does. For these reasons 
we would like to call it a "Bragg glass". Note that this is a much stronger property that 
the so-called "hexatic glass"E3EJ since hexatic order in the elastic limit is a straightforward 
consequence of the absence of dislocations. In the Bragg glass, two important length scales 
control the crossover towards the asymptotic decay, a consequence of the fact that the 
disorder varies at a much smaller scale than the lattice spacing a. When the mean square of 
the relative displacement {[u{x) — n(0)]^) of two vortices as a function of their separation x 
is shorter than the square of the Lindemann length = (m^), the thermal wandering of the 
lines averages enough over the random potential and the model becomes equivalent to the 
random force Larkin model. At low enough temperature. It is replaced by the correlation 
length of the random potential C,o, which is of the order of the superconducting coherence 
length. In that case the crossover length is When the relative displacement is larger 

than the correlation length of the random potential but smaller than the lattice spacing a, 
this is the random manifold regime where each line sees effectively an independent random 
potential. When the relative displacement becomes larger than the lattice spacing, one enters 
the asymptotic quasi-ordered regime. This occurs for separations of order ^. In general the 
two lengths Rc and ^ are widely different. The theory developed here can be applied to any 
elastic system in the elastic limit ^ ^ a. In relation with experimental systems we focus 
particularly on the triangular Abrikosov lattice in c? = 2 + 1, point vortices in thin films and 
magnetic bubbles in (i = 2 -|- 0, and will also mention lines in a plane d = 1 + 1. 

In this paper we will not treat topological defects quantitatively. Although a full descrip- 
tion of a lattice in presence of disorder should also include topological defects, their effect 
might not be as severe as commonly believed from misleading Imry-Ma type arguments. 
Indeed the fact that within the elastic theory quasi-long range order is preserved at large 
distances makes the system much more stable to dislocations. Since in ci = 3 the core energy 
of a dislocation increases as its size L, it is actually very possible that a phase without un- 
bound dislocations exists in = 3 in the presence of weak disorder. Indeed, Bitter decoration 
experimentsill at the highest fields available, about 70 G for these low-fields experiments, 
show remarkably large regions free of dislocations. In recent neutron experiments^^ it was 
shown that the degree of order depends on the way the system is prepared. A more perfect 
lattice with a smaller number of dislocations was prepared by first driving the system at a 
velocity high enough for translational order to heal, and then slowing it back down to zero 
velocity. It is thus conceivable that in d = 3 the presence of dislocations is overall a non- 
equilibrium feature. In two dimensions {d = 2 + 0), dislocations are energetically less costly 
and will probably appear at large scales, although this has not yet been firmly established. 
However, as we will discuss here, the length scale between unbound dislocations can be 
much larger than ^ in a low temperature regime. In that regime the main cause of the decay 
of translational order is elastic deformations due to disorder. 

The paper is organized as follows: For convenience we have separated the mostly technical 
sections ( |III| ) and (|V|) from the one discussing applications to physical systems ( [TVD and 
(|VTD. In Section (|ID, we introduce the model and derive the correct elastic limit. Simple 
dimensional arguments a la Fukuyama-Lee are given to identify the relevant length scales. In 



4 



section (|III|) we apply the Gaussian Variational Method to a simplified isotropic version of the 
model. Thus the method can be exposed without being obscured by unnecessary technical 
complications specific to real vortex lattices such as anisotropy and non-local elasticity, while 
the essential physics is retained. This section contains most of the technical details and 
methods used. We explain why a previous application of the variational method by BMY 
led to erroneous conclusions about the behavior at large scales. In Section (^V]) we apply the 
theory to ci = 2 + 1 solids such as the vortex lattices using a realistic elastic Hamiltonian. We 
compute the translational order correlation function Cxoir) with the full crossover between 
the three regimes in distance. We discuss the experimental signatures for decoration and 
neutron diffraction experiments. In particular the results of a comparison between decoration 
images and theoretical predictions are mentioned and detailed predictions are made for the 
neutron diffraction intensities. We then give a simple physical interpretation of the various 
regimes in distance and also argue that dislocations are less likely to appear as commonly 
believed. In Section (|Vp we apply the functional renormalization group in d = 4 — e, and 
compare its findings with those of the variational method. In Section (|V| ) we examine two- 
dimensional systems, for which thermal fiuctuations play a more important role. We first 
apply the variational method which shows that below a critical temperature Tc the system 
is glassy with logarithmic growth of displacements. This is compared to predictions of the 
renormalization group in d = 1 + 1. We then give a physical discussion of what should be 
expected for ci = 2 + systems, such as magnetic bubbles, where dislocations have to be 
considered. Conclusion can be found in Section |VI1| . Finally the bold and brave can look at 
the appendices where most of the most tedious technicalities are relegated. 
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II. DERIVATION OF THE MODEL AND PHYSICAL CONTENT 



A. A General elastic Hamiltonian 



In the absence of disorder the vortices form, at equihbrium, a perfect lattice of spacing a 
whose sites are labeled by an integer i and position will be denoted by Ri. Since we want to 
apply this theory also to a lattice of vortex lines, we consider the more general case where 
the Ri are m-component lattice vectors and there is in addition d — m transverse directions 
denoted by z so that the total spatial dimension is d. Throughout this paper we will denote 
the d dimensional coordinates by x = (r, z) and similarly the Fourier space (momentum) 
coordinates by g = {q±,qz)- For example, the Abrikosov lattice corresponds to m = 2 and 
d = 3 and z is along the direction of the magnetic field. The displacements relative to 
the equilibrium positions are denoted by the m component vector u{Ri,z) = Ui{z). For 
weak disorder (a/^ -C 1 where ^ is defined below ) and in the absence of dislocations, it is 
legitimate to assume that u{Ri, z) is slowly varying on the scale of the lattice and to use 
a continuum elastic energy, as a function of the continuous variable u{x). We consider the 
simple elastic Hamiltonian: 

^el = ^ E j^^ 7^^a(9)'^'a/3(g)M/3(-g) (2.1) 



where = 1, ..m labels the coordinates, and BZ denotes the Brillouin zone. The are 
the elastic matrix. Such an elastic description is valid as long as the relative displacement of 
two neighboring points remains small i.e. \u{Ri) — M(i?j+i)| <^ a, but does not suppose the 
individual displacements themselves to be small. We differ the study of the realistic elastic 
Hamiltonian ( |2.1| ) until sections |V| and [V|, and in order to illustrate the method in sections 
|I| and |ITl| we use the fully isotropic elastic Hamiltonian 



H.i = lj d'xiVuix))' (2.2) 



corresponding to the case ^aisiq) = cq^^ap- 

In the limit where many weak impurities act collectively on a vortex the disorder can 
be modeled by a gaussian random potential V{x) with correlations: V{x)V{x') = A(r — 
r')5{z — z') where A(r) is a short range function of range (the superconducting coherence 
length) and Fourier transform Ag^. The other limit, corresponding to a few strong pins, can 
be modeled by a Poissonian distribution and will not be considered here. Since the density 
of vortices at a given point is given by 

p{x)=J25{r-Ri-u{Ri,z)) (2.3) 

i 

the total Hamiltonian is therefore 

H = H,i + J d'^xV{x)p{x) (2.4) 
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The simplest way to recover translational invariance, is to use the well-known replica 
tricki. amounts to introduce n identical systems by replicating the original Hamil- 

tonian. It is then possible to average over disorder, the proper quenched average being 
recovered in the limit n ^ 0. After replicating, the interaction term becomes 




where a, b are the replica indices. (|2.5|) can be rewritten 

^Pin = E / d'^^d'^x' E A(a; - x'M^ - z'M^ -R^- <{z))5{x - Rj - u]{z)) 

a,b i,j 

= E E / rf'~™^A(i?, - i?, + u^z) - u%z)) (2.6) 

a,b i,j 

Where a,b = l,...,n are the replica indices. As we show in the following section it is 
extremely important to keep the discrete nature of the lattice in ( p.6|) , and the continuum 
limit of i^Din should be done with some care. 



B. Decomposition of the density 

Using the form ( |2.6| ) for the Hamiltonian in term of the displacement fields Ui is rather 
cumbersome. ( |2.6|) leads to a non local theory, even in the limit where the disorder is 
completely uncorrelated A{r — r') = S{r — r'). Indeed, vortices belonging to two different 
replica sets can be a priori at the same point in space r = Ri + u^{z) = Rj + m^(^) while 
having very different equilibrium positions Ri ^ Rj. This occurs when the displacements 
of the vortices are large enough. Since Ri, the equilibrium position of the vortices, have 
clearly no physical significance, except as an internal label, it is much more convenient to 
use instead a label that is a function of the actual position of the vortices. This can be 
achieved by introducing the slowly varying field 

(f){r, z) = r — u{(j){r, z), z) (2.7) 

Such a field will allow the continuum limit of (|2.6| ) to be taken easily. For each configuration 
of disorder, or alternatively for each replica set, one introduces a different field 0. Such a 
labeling is always exact when the transverse dimension m is m = 10. In more than one 
transverse dimension, this representation assumes the absence of dislocations in the system. 
In a self consistent manner, we will justify a posteriori both assumptions of elasticity and 
absence of dislocations using the solution in sections ( [111 L)| ) and ( [IV (J| ). 
Using (f){x) the density can be rewritten as (see appendix 0) 

p(x)=podet[9,0^]5:e*^-^(^) (2.8) 

K 

where the K are the vectors of the reciprocal lattice and po is the average density. In the 
elastic limit one can expand ([2.8[) to get 
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p(x) ~ po(l - + e'^'^pKix)) (2.9) 

where 

p^(x) = e-*^"^'^^^)) (2.10) 

is the usual translational order parameter defined in terms of the reciprocal lattice vectors 
K. Expression ( |2.9| ) respects the periodicity of the lattice i.e. is obviously invariant by a 
global translation u ^ u + a. Another advantage of the decomposition (|2.9| ) is that the 
various Fourier components of the density relative to the periodicity of the ordered lattice 
appear clearly in H^i^ 

/ dxV{x)p{x) = — po / dxV{x)daUa + Po dx Y^ V-Kix)pKix) (2-11) 

where 

Vk{x) = V{x)e-'^' (2.12) 

is the part of the random potential with Fourier components close to K. Since the energy 
is invariant when changing u ^ m + a, u itself cannot appear in the expression ( |2.11| ), but 



only 9m, and in principle higher derivatives, or a periodic function of u are allowed. 

The first term in the right hand side of ( p.ll| ) is the part of the deformation of the lattice 
that couples to the long wavelength of the disorder potential. It results in an increase or 
decrease of the average density in regions where the potential is favorable or unfavorable. 
The second term couples to the higher Fourier components of the disorder. The average 
density is not affected but the lattice can be shifted so that the lines sit in the minimum 
of the disorder potential. In the usual elasticity theory, one takes the continuum limit for 
the displacement field and assumes that the density itself is smooth on the scale of the 
lattice. This allows to keep only the gradient term in ( p.ll| ). Here although it is possible to 
take the continuum limit for the displacements u since they vary slowly on the scale of the 
vortex lattice ( Vm -C 1) , it is imperative to retain the discrete nature of the density. This 
is because the scale at which the disorder varies (for superconductors it is comparable to 
the scale of the real atomic crystal) is usually shorter than the lattice spacing of the vortex 
lattice itself. 

If one uses the representation ( |2.9| ) of the density and ( ^.6] ) and discard spatial averages 
of rapidly oscillating terms, the replicated Hamiltonian becomes, in the isotropic case: 

H,s = y d''xiVuix)f (2.13) 

•' a,b K^O 

To be rigorous the last terms in (|2.13|) should be written in term of u{(j){x)) rather than 
u{x), but this has no effect on our results. It leads only to corrections of higher order in Vu 
which we neglect since we work in the elastic limit a/^ -C 1. The Hamiltonian (|2.13| ) will 



be our starting model, and from now on we absorb the coefficient pg in A^, PqAk —>■ Ak- 
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A general property of the Hamiltonian ( p.l3| ) is the invariance of the disorder term 
under the transformation Ua{x) — >• Ua{x) + w{x) where w{x) is an arbitrary function of x. 
This statistical invariance guarantees that the elastic term in ( p.l3| ) is unrenormalized by 



disorder. Note that in the original non local model (|2.6|) this symmetry is only approximate, 
and indeed one would find there a small (of order (a/^)^ ) and unimportant renormalization 
of the elastic coefficients by disorder. 

The principal quantities of interest are the mean squared relative displacement B{x) of 
two vortices, averaged over disorder, which is determined by the correlation of u diagonal in 
replicas 



B{x) = -{{u{x) - u{m = 2tJ ^(1 - cos(gx))G(g) (2.14) 
and the translational order correlation function Ck{x) 



Ck{x) = {pMpKm (2.15) 

In the gaussian theory that will be considered below, the two correlation functions are simply 
related by 

Ck{x) = e-^""^'^ (2.16) 



The Hamiltonian ( p.l3[ ) can be applied directly to study quantum models with a time 



dependent disorder. A more physical disorder for quantum systems would be only space 
dependent. This would correspond to correlated disorder in one (the "time") direction for 
the classical system and can be studied by the same method than the one used in this 
paperE3. For completeness we also give here in appendix |^ the connection between the 
quantum mechanics of interacting bosons and fermions and an elastic system in c? = 1 + 1 
dimensions. 



C. Dimensional Arguments 

Before starting the full calculation, let us estimate the effects of the different terms in 
( [^.11| ) in a way similar to Ref. ^ In the presence of many weak pins, u cannot distort 
to take advantage of each of them, due to the cost in elastic energy. One can assume that 
u varies of ~ a over a length ^ ^ a. The density of kinetic energy is ~ c(a/^)^, where c 
is an elastic constant. The various Fourier components of the disorder will give different 
contributions. The long wavelength part of the disorder gives 

^Ji'o ~ Po/ d''xVix)d^Uaix) ~ Al^^/e-""/^ (2.17) 

For the higher Fourier components the disorder term 

H^}^j^ = PoJ d'^xV^(x)e*^^e-*^"(^) (2.18) 
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can be estimated over the volume ^'^ as e~^^^ df'xV {x)e^^^ . This sum can be viewed as 
a random walk in the complex planed and the value of u adjusts itself to match the phase 
of the random potential. Therefore the gain in energy density due to the disorder term is of 
order 

/e'/' (2.19) 

Optimizing the gain in potential energy gain versus the cost in elastic energy determines 
^. One can therefore associate to each Fourier component a length scale above which the 
corresponding disorder will be relevant and destroy the perfect lattice 

Cg^o ~ a (c^a7Ao) if < 2 (2.20) 

e,.^~a(cV/A^)'/^'"'^ (2.21) 

The g ~ component of the disorder is relevant only for d < 2 and the second term in 
( ^.l^j| ) can be dropped for d > 2 if one is interested in the asymptotic regime. In fact, the 



g ~ part of the disorder can be eliminated exactly from ( |2.13|) and leads only to trivial 
redefinitions of the correlation functions. One can perform a translation of the longitudinal 
displacement field u by 

Ua{x) Ua{x) + /„(x) (2.22) 

where the Fourier transform of / is 

= (2.23) 



The translation ( |2.22D when performed on the replicated form ( ^.13| ) eliminates the long 



wavelength term but does leave the cosine term invariant since it is a local transformation. 
Note that such a transformation is only possible due to the fact that the various Fourier 
components of the disorder are uncorrelated. The mean squared relative displacement B{x) 
becomes 

As expected the additional term produces only a subdominant finite correction above two 
dimensions. In the following we simply set Aq = 0. 

As is obvious from ( |2.20| ), higher Fourier components Vqr^x disorder the lattice below 
d = 4. We will now examine the effect of these Fourier components more quantitatively. 
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III. VARIATIONAL METHOD 



We now study the Hamiltonian ( p.l3| ) using the variational method introduced by Mezard 
and Parisiii. Hamiltonians with more realistic elastic energy terms, directly relevant for 
experimental systems, will be considered in sections ( fVl) and (|VlD. 



A. Derivation of the saddle point equations 

We now look for the best trial Gaussian Hamiltonian Ho in replica space which approx- 
imates ( p.l3| ). It has the general formil: 



^0 = ^1 ■^^Gj{q)Ua{q) ■ Ubi-q) (3.1) 

where the [G~'^]ab{q) is a n by n matrix of variational parameters. Without loss of generality, 
the matrix G~b (q) can be chosen of the form G~b = cq^Sah — o^ah where the self energy aab 
is simply a matrix of constants. The connected part is defined as G~^{q) = Y.bG~h{q)- We 
obtain by minimization of the variational free energy F^ar = -^o + {Heff — Hq)hq the saddle 
point equations: 

^...-E^A-=^-*-<-. (3.2) 

where m is the number of components of u. One defines 

Babix) = ^{[uaix) - UbiO)]^) (3.3) 



T I -^AGaaiq) + GM -2cos{qx)GM) 
J (27r) 



Note that the connected part is unchanged by disorder, a direct consequence of the statistical 
symmetry of (|2.13| ) noted above. Two general classes of solutions can exist for (|3.2|) . One 



preserves the symmetry of permutations of replica, and amounts to mimic the distribution 
( thermal and over disorder) of each displacement mode u{q) by a simple Gaussian. The 
other class, which is a better approximation in the glassy phase, breaks replica symmetry and 
approximates the distribution of displacements by a hierarchical superposition of gaussians 
centered at different points in spacec3. Each gaussian at the lowest level of the hierarchy 
corresponds to a different metastable "pinned" position of the manifold. 



B. Replica Symetric solution 

Let us first examine the replica symmetric solution Ga^bio) = G{q) and Ba^b{x) = B{x) 



Using ( |3.3| ) one has 

B{x = 0) = 2t| ^,GM (3.4) 
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For d < 2, B{x = 0) is infinite, and the off-diagonal part of aab is zero. The K ^ Fourier 
components of the disorder do not contribute. This solution turns out to be the correct 
solution for c? < 2, as shown in appendix |C[ This can be explained physically by the fact 
that for d < 2, thermal fluctuations are strong enough to disorder the system. 
For d > 2, G is given by 

'5(')-^^ + W^SA.A-V-4/^ (3.5) 

It is the Lindemann length and measure the strength of thermal fluctuations. It is defined 
as: 



where 1/5"^ = 2'^'^^n'^^'^T[d/2]. Due to the term in in ( ^.51 ), the relative displacement 
correlation function grows as 

B{x) ~ x^-'' (3.7) 

The replica symmetric solution is therefore equivalent at large distances to the Larkin result 
based on a model of independent random forces acting on each vortex. As explained in the 
introduction this solution does not contain the right physics to describe the long distance 
behavior. In this variational approach this shows by the fact that the replica symmetric 
solution is unstable towards replica symmetry breaking for 2 < d < A. This can be checked 
from the eigenvalue A of the replicon mode@. 



^_j_^A^^-.-../j-c.,„ Ij- (3.8) 



K 



m J (2ti 



A negative eigenvalue A indicates an instability of the replica symmetric solution. We 
introduce a small regularizing mass in Gc- Gc{q)~^ = cq^ + /i^ and take the limit fi —>■ 0. 
It is easy to see from (|3.8|) that for d < 2 the replica symmetric solution is always stable 
(see also appendix |y). In that case disorder is in fact irrelevant, due to the strong thermal 

fluctuations. For d = 2 the condition becomes < 1 for small fi. Thus there 

is a transition at T = Tc = Attc/Kq between a replica symmetric stable high temperature 
phase where disorder is irrelevant and a low temperature (glassy) phase where the symmetric 
saddle point is unstable. We will examine the physics in = 2 in details in section ^I[ For 
2 < d < A the replica symmetric solution is always unstable and disorder is therefore always 
relevant. 



C. Replica Symmetry breaking for 2 < d < 4 

Since for 2 < (i < 4, the replica symmetric solution is unstable, to obtain the correct 
physics one has to look for a replica symmetry broken solution. We will focus here on 
the case 2 < ci < 4, the d = 2 case being discussed in section |VT[ Following^ we denote 
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G{q) = Gaail), similarly B{x) = Baa{x), and parametrize Gabil) by G{q, v) where < v < 1, 
and Bab{x) by B{x,v). Physically, v parametrises pairs of low lying states, in the hierarchy 
of states, as described inlHj v = corresponding to states further apart. The saddle point 
equations become: 

^(^)= V^i^V^^'^(°''') (3.9) 
^ mT 

where 

i?(0, ^) = 2T / ^,{G{q) - G{q, v)) (3.10) 

B{0,v) corresponds physically to the mean squared relative displacement between the posi- 
tion of the same vortex (x = 0) when the manifold is in two different low lying metastable 
states. The large distance behaviour of disorder-averaged correlators will be determined by 
the small v behaviour of B{0,v), 

As we will show in section |111 C 2 , to discuss the large distance behavior x ^ ^ it is enough 



to keep the smallest reciprocal lattice vectors with = Kq in ( |3.9| ) since i?(0,f) ^ a?. 
We will thus first study a single cosine model obtained by keeping only K = Kq 

Hcs=y d'^xY^iVu'^ix))'' - ^ cos(iro(w"(x) - u\x))) (3.11) 

Each time we consider this particular model, e.g in the following subsection ( [III C 1|) , we will 
denote by A = NbAxo where A*";, is the coordination number, i.e the number of vectors Kq 
with minimal norm. 



1. Asymptotic behavior (single cosine model) 

We look for a solution such that cr{v) is constant for v > Vc, fc itself beeing a variational 
parameter, and has an arbitrary functional form below Vc- The algebraic rules for inversion 
of hierarchical matrices!! give: 

^ ' ^ ' Jv J {27rr{G,{q)-^ + [a]{w)f ^ ^ 

where [o']{v) = ua{v) — Jq dwa{w) and 

B{0, v^) = / \ (3.13) 

In that case, taking the derivative of ( pl9|) ( keeping only K=Kq ) with respect to v, using 
[cr]'(f) = va'{v), ( |3.12|) , and (|3.9| ) again one finds 
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Since the integral is ultraviolet convergent, we have taken the short- distance momentum 
cutoff A = 27r/a to infinity, a limit discussed below. Derivating one more time one gets for 
the effective self energy: 

[cr]{v) = {v/v^fl' (3.15) 
where Vq = 2K^TcdC-'^^'^ /{A - d) and 

~ i (2^ l^^^TTy ~ 2'^+isin(rf7r/2)r(t^/2) ^^'^^^ 

with Crf=3 = l/(87r), Cd=2 = l/(47r). 

The behavior of [(y]{v) controls the scaling of the energy fiuctuationsta AF oc oc T/v, 
with the scale L, and the large scale behavior is controlled by small v. (|3.15|) thus gives an 
energy fluctuation exponent 6 = d — 2. 

Using ( |3.15|) in (|2.14| ) one can now compute the correlation functions. Larger distances 
will correspond to less massive modes, and one obtains 

{{u{x) - u{0))^) = 2mT I ^(1 - cos(gx))G(q) (3.17) 



GM = ij(l + /' ^ ~ (3,18) 



with Zrf = (4 - dY/n'K^Sd) and 1/^^ = 2'^'^n'^/^T[d/2]. Thus for 2 < < 4 we find 
logarithmic growthEaS: 



2m 



-m(0))2) = A,log|x| (3.19) 

with Ad = 4: — d. Note that the amplitude is independent of temperature and disorder. 

The solution ( ^.15 ), is a priori valid up to a breakpoint Vc, above which [a] is constant, 
since cr'{v) = is also a solution of the variational equations. To obtain the behavior at 
shorter scales for the single cosine model, we need to determine the breakpoint Vc- For v > Vc 
using [cr](f) = S one can rewrite ( p.l5| ) as 

M(^) = S(^)i (3.20) 

with Vc = VoTi~ . Using (|3.9|) and ( p.l4| ) the equation determining S is: 

= ^e-'^^oBio..) (3.21) 

in terms of the nonuniversal quantity B{0,Vc)- 

^ f d'^q 1 
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One can define a length / such that c/~^ = S. Since x <^ / corresponds to v Vc where 
is a constant and the solution is similar to a replica symmetric one. / is therefore the 
length below which the Larkin regime will be valid. When / 3> a, one finds B{0,Vc) ~ ^r^- 
For instance in = 3, B{0,Vc) = /t^(1 ~ (a/(27r/)) arctan(2'7r//a)). Equation ( p.21|) can be 
rewritten 

which is equivalent to Arepiicon(o") = 0. Assuming / ^ a one finds 
with 

e = (mcVAiro^Q)^/^^-'^) (3.25) 

and 



2irgrc^ 

(4 - d)cl^~^ 



(3.26) 



Note that although the breakpoint Vc 0, when T — > 0, the length I which is associated to 
the transition between the two regimes remains finite. 

For the single cosine model, the characteristic length / below which the replica symmetric 
part of the solution ([o"]('u) = S) determines the physics, is equal to ^ the length for which the 
relative displacements are of order a. For this model one has a direct crossover between the 
Larkin regime and the logarithmic growth of the displacements. This is to be expected since 
the disorder is here characterized by a single harmonic. It has therefore no fine structure for 
distances smaller than a, the lattice spacing. This will not be the case any more if higher 
harmonics are included. The disorder will be able to vary strongly for distances smaller than 
a, and one expects / and ^ to be different, and a third regime to appear in between: the so 
called random manifold regime. 



2. Study of the crossover 



We study now the full model (|2.13|) . Since this model contains all the harmonics of the 



disorder, it can describe correctly the short distance regimes. In particular, we will examine 
here the crossover from the random manifold regime to the logarithmic one. 

In order to rewrite the equations in term of dimensionless quantities, we introduce the 
rescaling 

a(v) = -^s{v/v^) (3.27) 
,2 
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as will be obvious later ^ is the crossover length between the random manifold regime and 
the logarithmic one, and corresponds to the value of v for which the crossover occurs. One 
chooses ^ and such that ( |3.9| ) and ( p.l2| ) become in terms of the dimensionless quantities 



b{y)=b{y,)+ Tdy- "'^^^ 



(3.28) 
(3.29) 



where yc = f c/ and the integration over momentum in ( |3.12|) has been performed. We have 
introduced the dimensionless variable p such that K = 27Tp/a. When using the definition 
( p.28| ) one gets 



4 2 

ma c 



27r2 f2Tcdfl^\ (a 
«^ I c / U. 



(3.30) 



d-2 



It 



a 



d-2 



Thus is always very small compared to 1. In ( p. 301) , for simplicity, we have assumed that 
all Ak have the same value A = Akq- 

The equations (|3.28| ) and ( |3.29|) can be solved in a parametric form. We introduce the 
variable z = b{y) and define 



Hiz) 



(3.31) 



It is possible to keep different Ak, for instance ~ exp{—K'^^l) to describe the effect of 
the finite correlation length of the random potential, by just modifying the functions H and 
h = -H' to: 



(3.32) 



Using the variable z and taking the derivative of ( p.29|) , the equations ( |3.28|) and ( p.29|) 
become 



s{y) = H{z) 

h\y) - 



s'{y) 



[s\{y)i*-d)/2 

Taking the derivative of ( p.33|) one gets 

s{y) = Hiz) 
[sjiy)"^ = h{z) 



(3.33) 
(3.34) 



(3.35) 
(3.36) 



Finally, using y 



= ^ ^6 obtain the solution in a parametric form 
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4-d 



[s]{y)^ = h{z) (3.37) 

2 2d 6 

h\z)h{z)'-^ (3.38) 



with Zc = h{yc) < z < oo. 

Let us examine first the various asymptotic behaviors of the solution ( |3.37| ). As will be 
obvious later, large z correspond to large scales and small z to small scales. At large z only 
the smallest p contributes in the sum (|3.31D for h{z) giving 



h{z) = 2me ^ for a square lattice (3.39) 
h{z) = ^e~^^^^ for a triangular lattice 



In that case the high harmonics are irrelevant and ( p.37 - 3738| ) give back formula ( p.l5 ), for 



the single cosine model. One recovers the quasi ordered large distance logarithmic regime. 

We now study the behavior at small z. In that case all harmonics must be kept and it 
is convenient to use the following duality transformation of formula (|3.31|) 



/W=Ee-'' = ^i: (3.40) 

where the vector p have been defined above. Q is the volume of the unit cell in the space of 
the vector p [Vt = 2/ for a triangular lattice and 1 for a square lattice). The vector R are 
the reciprocal vectors of the p which themselves are normalized in units of 27r/a. The R thus 
correspond to the original lattice with a spacing unity. For small only i? = contributes 
and 

I{z) ~ - ( - ) (3.41) 

therefore 



Vt \z 



which gives 



, , , 7r™/2 m m + 21 , , 

hiz) ~ — ^ r 3.42 



[s](2/) = Sy2/^«M (3.43) 



where the fluctuation energy exponent of the random manifold regime is ^ = {2d — 2m + 
dm) /(4 + m) and the amplitude 



~ll 4 j V4 + mj 



(3.44) 



The equations ( [^.14D and ( |3.18| ) once rescaled using ( ^.3U| ) give 



2 



B{x) = —h{-) (3.45) 
2-K^ 4 
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and 



with the useful intermediate formula: 

dy y"" n/a 







y'^ 1 + y^ sin vr/a 
Using the asymptotic expression 

d'^q 1 



(3.47) 



{2TTyq 



1 - cos(gx)) = Id.yX^" (3.48) 



with 



TT 



2'^+2^r(d/2 + z/)r(l + z/) sin(z/7r) 



(3.49) 



one gets 



Kx) ~ ( . ... s^'--^^/.,. I (3.50) 

VQ2sm(7r6'RM/2) / 

with 2z/ = 2 + — c^- Thus the exponent entering in the relative diplacement growths is 
z/ = with z/ = l/6form = 2. This corresponds to the random manifold regimeSil. 
In this regime each vortex is held by the elastic forces of the other vortices and sees an 
independent random potential. The mean squared displacement grows more slowly than in 
the Larkin regime (the exponent is 1/3, compared to 1 for the Larkin regime in c? = 3). In 
c? = 3 and m = 2 one gets ^rm = 4/3 and for a triangular lattice in c? = 3 

the amplitude in ( ^.5(J| ) is ~ 2.3817. 

As for the single cosine model of section [III C 1| , the solution (|3.43|) is valid up to a 



breakpoint Vc above which the self energy [o']{v) is constant. This corresponds to scales such 
that B{x) is smaller than and ^q. One then recovers the replica symmetric propagator 
G{q) ~ 1/g^ for q^ ^ [cr](fc), and Larkin's model behavior. To compute the crossover and 
to determine the breakpoint Vc = yc'v^^, we proceed similarly to section [111 C 1| . The equation 
determining S is now 

where 5(0, fc) is given in ( p.22[ ). Note that keeping the correlation length of the disorder 
using Ax = A exp(— 1/2X^^0^) amounts to change 5(0, Wc) into 5(0, fc) + and thus /f, 
into + ^0^. We will thus take = A keeping in mind this change. 
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Solving (p.52|) and using the small z expansion of h(^z) one gets for the lengthscale / such 
that S = c/~^ (assuming / ^ a and It -C a) 




where 



(3.53) 



Cs = {2n-y/^^^\ , ^J V(4-.) (3.54) 



^ W2m(m + 2)' 



The breakpoint Vc can obtained using (|3.371 ) with the value z = Zc = 27r^/|,/a^. This gives 
using ( p. 30] ) 




(3.55) 



•.4^)"^f?r^^ (3.56) 



which leads to 



The characteristic length / separates the Larkin regime from the random manifold regime, 
and is in that case much smaller than ^. Lowering the temperature reduces the range 
over which the Larkin regime occurs. This is because the thermal wandering responsible 
for smoothing the random potential on a scale It decreases. The relative displacements of 
two vortices separated by / is of order Ma.x{lT,C,o), giving B{1) ~ Max(/T, ■Co)^- As T — >• 0, Z 
becomes identical to the Larkin-Ovchinikov length Rc. Using the expression of B{x) ~ Cx^*" 
in the random manifold regime and the additional relation B{^) ~ one recovers the 
expression ( p.53| ) for /. When f c = 1 the Larkin regime disappears. This occurs when / ~ a. 
The criterion / ^ a for which the Larkin regime exists is equivalent to 

(3.57) 

and corresponds therefore to extremely weak disorder and intermediate temperatures. The 
absence of a Larkin regime means that the disorder-induced relative displacement of two 
neighbors in the lattice is already larger than Max(/7-,^o)- 




3. Crossover in d = 3 

In d = 3, it is possible to solve the equations describing the crossover analytically, and 
thus to obtain the full crossover function between the random manifold and the quasi ordered 
regime. We will examine d = 3 and m = 2. for the model (|2.13|) . Such a case is physically 
relevant for the case of vortex lattices. The crossover length are given by ( |3.30| ) 
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Using p.45|) one has 



b{x) = - 

C3 



(2^ 



3A 



(3.58) 



'1 — cos(A;x)) 



V'^ dy 



[s]iy) 



k^{k^ + [s]{y)) 



Performing the angular integration over momentum in (|3.59| ) we find 



b{x) 



2tt'^cs Jo 



dk(l — —— sm(k\x\)) 
k\x\ 



dy [s]iy) 



(3.59) 



(3.60) 



where we have extended the integral over y to infinity, assuming -C Vc or equivalently 
It CL, in which case there is a wide random manifold regime. Performing the remaining 
integration over k one gets 



b{x) 



Arccs JO ?/ \ \x\ 



-\x\[s]'^Hyh 



(3.61) 



Using the parametric solution (|3.37|) and ( |3.38|) for [s] {y) we obtain the final expression 



b{x) 



°^ h"{t)h{t) 



dt- 



h'{ty 
1 , 



f{xh{t)) 



X 



(3.62) 
(3.63) 



Expression ( p.62| ) gives the full relative displacements correlation function as a function of 
distance. To recover the asymptotic expression of section |III C 1| , for large distance x, one 
notices that in ( |3.62| ), h{t) ~ y4e~"* as shown in ( |3.39D for large t. Thus the large x behavior 



(3.64) 



will be controlled by small t. One obtains the asymptotic expression 



= - 

where we have used [h/h' 



/(oo)(log(Ax) + 1) + / dz\og{l/z)f\z 



oo 




b{x 



1/a. Using ( p.62| ), one finds 
1 



a 



[\og{Ax) + 7] 



(3.65) 



where a = [KQa/{2TT)]'^ and 7 ~ 0.57721 is the Euler constant. This implies that the 
translational order correlation function 



{e 



iK(i-u(x) —iK(s-u{Qi) 



behaves for large 



A\x\ 



(3.66) 



(3.67) 



One recovers the power law behavior of section |111 C 1| as well as the amplitude. The in- 
termediate distance behavior will be examined in more details for more realistic elastic 
Hamiltonians in connection with vortex lattices, in section 0. 
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D. Self-consistence of the physical asumptions 



Finally, for our solution to be valid, one has to check self-consistently that even in the 
presence of disorder, the basic assumption that elastic theory was applicable remains valid. 
One has therefore to check that ^((Vn)^) <^ 1. 

For simplicity we will make the analysis for the single cosine model (for which / ~ C,) 
but similar results can be derived for the full Hamiltonian. Using the variational solution of 
section 1111 (J 11 one obtains 



m 



2T 

c 



cg^ + S(t>/t>c)^/^ Vc 



1) 



cq^ + S 



(3.68) 



Using (|3.6| ) one gets, for the case / ^ a 



m 



■{(yuf) 



/2 



cVr A — d 



+ 



, 2ttIt s2d — 2 

'a d 



(3.69) 



where the last contribution is due to thermal fluctuations only. Replacing Vc by (|3.26| ) one 
finds 



((V«)2) 



;sin(7r(c/- 2)/2) / a V''^ ^^t^^t^ 



7r{d-2 



\2-kIJ 



d-2 
d 



When d 2 one cannot neglect S in the denominator of ( 3.68 ), and the expression (|3.70 ) 
becomes 



log(//a) 



(3.70) 



As is obvious from ( p.70|) , one has always ((V'u)^) <^ 1, provided that / 3> a or equivalently 
using ( |3.24D for the single cosine model, provided that one is far from the melting temperature 
It ^ CL and that the disorder is weak ^ 3> a. In that case one can indeed use an elastic 
theory in the absence of dislocation, even in the presence of disorder, and our solution is 
valid in such a regime. 



E. Comparison with BMY 

The previous application of the variational method by BMY§'i led to the erroneous 
conclusion that the fluctuations are enhanced at large scale. They find for d = 3 B{x) ~ 
x^f^ instead of logarithmic behavior found here. Although they want to describe the same 
physical situation as the one studied here, they in fact consider a model which turns out to 
be fundamentally different, in which each vortex sees a different disorder. In their model 
the random potential is also dependent on the line index i such that 

V{Ri,r,z)V{Rj,r',z') = A6{r - r')6{z - z')\Ri - Rj\-^ (3.71) 

This amounts to introduce an extra disorder in the original model (|2.4| ) with correlations 
decaying as l/\Ri — Rj\^. Then BMY retain only the long wavelength part of this disorder, 



21 



which indeed for a fixed A > dominates the contribution of higher harmonics. They 
then look at the hmit of the exponents when A — >■ 0. The result they obtain with this 
procedure is incorrect (although their derivation is technically sound) and comes from the 
following artifact: by assuming that different lines see different random potentials, they 
make it possible to optimize the pinning energy by a global translation of the whole lattice. 
In that case the pinning energy will obviously be dependent on u, even for a uniform u. On 
the other hand for the genuine disorder, which is only dependent on the space position of 
the lines, it is obvious that a translation by one of the vectors of the lattice cannot change 
the energy, and therefore the g ~ part of Hp[^ cannot depend on u but only in du. By 
regularizing the integrals with a disorder dependent on the line index they introduce an 
extra and non-physical disorder which is relevant and changes the long-range behavior of 
the correlation function compared to the physical case. Indeed there is a crossover length 
associated with this disorder above which the long distance behavior is the one given by 
BMY. Below this distance the vortices all see the same random potential. To recover the 
physical model one has to take A — > and in that case C,x oo. 

In more mathematical terms, the variational method gives three types of contributions 
for the self energy as shown in appendix 

a{q, v) ~ cig2 + cae'^ + XB{0, v)-^-^/^ (3.72) 

The first term is the long wavelength contribution of the genuine disorder which is irrelevant. 
The second one is the higher harmonic contribution which is responsible for the logarithmic 
growth at large distances, and the third term is the long wavelength contribution of this 
extra disorder giving a B{x) ~ x^^^ for x > ^x. Only the third term was kept by BMY, 
artificially taking the limit A ^ in the exponent only but not in the amplitude of such a 
term. Note that if one takes the limit A — (which corresponds to the physical situation) 
before taking the limit x — oo one recovers that the g ~ part of the disorder does not 
play any role, and the amplitude they obtain vanishes. 

A simple Flory argument can be made to estimate the effect of the long wavelength part 
of the disorder on the displacements. This confirm that it is irrelevant above d = 2 (see also 
section |1I C|) . Let u be the typical relative displacement over a lengthscale L. The elastic 
energy cost is u^L'^''^ while the typical energy gain due to the disorder is 

^L'^ + uL'^-^ - oc (3 73) 

which comes from the change of density of the vortices. Since the vortex in the center are 
unaffected the gain of energy can come only from boundary terms. Balancing the two terms 
one finds u ~ l}-'^-'^)!'^ which is obviously irrelevant above two dimensions. 

In fact one can simplify the the saddle point equations oS@ by noting that the x 
dependence of B{x^ u) in these equations is unimportant, up to higher order terms in Vm. 
Such a calculation is performed in the appendix 0. One then recovers the local model 
( p.l3| ) which is simple and physically transparent enough to allow for the exact solution of 
section [111 C| . 
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IV. FLUX LATTICES 



A. Model 



The theory developed in section |T|, when specialized to m = 2 and d = 3, can be applied 
to describe the effects of weak disorder on the Abrikosov phase of type II superconductors. 
High Tc superconductors can be modeled by stacks of coupled planes. The system is therefore 
described by layers of two dimensional triangular lattices of vortices. The displacements u 
are two dimensional vectors, hence m = 2 (the vortex can only move within the plane). We 
denote by Ri the equilibrium position of the vortex labeled by an integer i, in the x — y 
plane, and by u{Ri, z) their in-plane displacements, z denotes the coordinate perpendicular 
to the planes and along the magnetic field. The total energy is: 

H = ^ J d^rdz[{cn - C6G){daUaf + CeeidaUfsf + Cu{d,UaY] + J d^rdzV{r, z)p{r, z) (4.1) 



where a, (3 denote in-plane coordinates. The Hamiltonian ( |4.1| ) is identical to ( p^.l| ) with 



^aM = G-''^{q)P^, + G-^^\q)Pt, where 



G-^^^{q) = c,,ql + ceegi (4.2) 
G-^^^{q) = cuql + cnql 



and 



PapiQ) = ^af3 - qo^qpl q± (4.3) 
PafsiQ) = Qaqp/ql 

are the transverse and longitudinal propagator. q±_ denotes the in plane vector, whereas qz 
is the out of plane component. Equation ( [4.1|) corresponds to a local elastic theory, but non 
local elasticity can also be considered at the expense of introducing q dependent coefficients 
c. This point will be considered in greater detail below. For the moment we restrict ourselves 
to dispersionless elastic constants. 

Weak point-like disorder such as oxygen vacancies, or defects introduced artificially in a 
controlled way, e.g by electron irradiationO, can be modeled by a gaussian random potential 
of correlation length of order ,^o- Here the disorder will be taken as completely uncorrelated 
from plane to plane A(x — x') = A(r — r')6{z — z'). Such a description will be valid as long 
as each pinning center is weak enoughso that the pinning length / (also called Rc in formula 
(51) of the Larkin-Ovchinikov papei^^) is much larger than the average distance between 
impurities. 

The disorder term in ( |4.1| ) is transformed into a form similar to ( ^.13| ). We can now use 
the methods of section |ITl| with the realistic elastic Hamiltonian ( |4.1|) to get the physical 
properties of a vortex lattice. Most of the theoretical calculation are confined in section [iV B| , 
whereas a simple physical interpretation of the results is given in section |IV C[ The experi- 



mental consequences are discussed in details in section IV D 
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B. Theoretical predictions 



One can then perform a variational ansatz identical to (|3.1|) but for the introduction of 



the longitudinal and transverse part G 



and 



Gai3,abiq) - G'^i,{q)P^p + G^^{q)P^p 
The saddle point equation ( |3.9D for the self energy now becomes 

A 



niT 



(4.4) 



(4.5) 



(4.6) 



K 



The correlation function B^p is defined by ( |3.3| ) with the replacement of G by Gap- 

Since B^s{x = 0, v) is a purely local quantity it is isotropic and B^s{x = 0,v) = 5^sB{x = 
0,v). This implies that aa/3 = o'{v)Sap, i.e. an isotropic self energy. Thus B{v) = ^[B^{v) + 
B'^{v)], where by definition B^''^{v) satisfy equation (|3.12| ) with respect to G^'^ . Integration 
over q leads to 



B(v) 



1,1 1 
T72(— + — . 

C44 Cu 



Bivc)+Tcd^i— + —) 



'{v)dv 



(4.7) 



The solution of (46) can trivially be obtained from the isotropic solution with the replace- 
ment 



C — > C44 

, 1 . C44 C44 , 

Q ^ Crf = Q-( + ^ 

2 C66 Cii 



(4.8) 



One can now compute the correlation functions 



Bal3{,x) = {[Ua{x) - M„ (0)] [M/3 (x) - M^(0)]) 

B^P:^^ir)+B^P^p{r) 



(4.9) 



where the longitudinal and transverse propagators have been defined in ( [4.3| ) and similarly 

to ( gTzD 



5^(r,2) 



2T 



(glf)2(l - cos(g_L'^ + qzz))G^{q) + 



(27r)3 

[1 - (glr)2](l - cos(gxr + g,z))G^(g) 



(4.10) 



and a similar equation for B^ obtained from ( [4.10| ) by permuting L, T. G'^'^ are defined 
similarly to equation ( |3.18| ) with the replacement of cq^ by G~^'^^ defined in ( [4.2| ). One 
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then rescales q and r to obtain isotropic integrals over momenta. Equation ( |4.10| ) then takes 
the form 



a C44 /C44 



27r2 'cii 



C44 f^T / /'^44 



z) + —F' [r 

Cll C66 




^)] 



(4.11) 



Those integrals contain the self energy [cr](f ) which is determined itself from equation 
4.7). These equations are rescaled similarly to ( p. 27 ) with the replacement (|4.8|) . This 
defines two crossover length ^ and for in plane and z directions, given by 



^C66 Cll ^ 



2fl4 1/2 3/2 
za C44 Cgg 

C44'' 7r3A(l + £^) 




(4.12) 



where, as we recall A is the disorder strength poV{q)V{—q) = A. 
Rescaling by the length ^ and ^2 one gets 



B'^ir^z) 



a ^C44 Cqq z C44 -j, r z 



+ 



27r2 Cll ^ V ^11 ^66 ^ ' 

where the functions F are given, for 2; = by 



(4.13) 



F'(r) 



1 



(27r)34 Jo y 



dy. 



[s]{y) / cos^(^)(l — cos(g_Lr cos(^))) 



1 



1 



%s]{y) I £q^dqz{l - cos\e)){l - cos(gxr cos(^))) 



qI + qI + + M(y) 
1 1 



(27r)3c^7o y^ 
performing the q^^qz integrations, one gets 



T± + qlqi + ql + [s\{y) 



(4.14) 
(4.15) 



1 dy 



47rc3 JO y"^ 



(4.16) 



where 



2 [s]i/2|a;|2 [s]V2|a;|2^ 
/r(x) = fi{x) - fL{x) 



s]l/2|; 



(4.17) 



a; a; 



Expressing again the [s] in term of the functions h, one gets 



. 2c 



11 IT,L 



27r2 Cll + C66 ^ Cll + C66 

where 6^'^ have an expression similar to (|3.62| ) 

dt 



2C66 ]jL,Ti^ / C66 ^ 



(4.18) 



/l'(t)^ 



(4.19) 
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with /^'"^ given by 



1 1 

- ^ + 



1 1 



2 X x^' 



(4.20) 



( |4.1S - O0| ) give the complete expression of the displacement correlation function for equal 
function of the distance in the transverse plane. 
In the large distance regime one obtains an expression similar to ( p.64|) , with the / 
replaced by f'^'^. This gives 



6^(x) 



1 

2a 
1 

2a 



log(y4x) + 7 + 
log(y4x) + 7 - 



(4.21) 



where a = [Koa / {2it)]'^ , 7 ~ 0.57721 is the Euler constant, and A = 32/3 for the triangular 
lattice. This leads for the B^'^ functions at large distance 



^0 



log{Ar/0 + 7 + 



e-^'"^ Cii — Cge 



C66 



■log(— ) 

2 Cii + Cqq 2(cii + Cqq) Cii 



(4.22) 



where = -1, = +1. Note that B = {B^ + B^)/2. 

It is interesting to note that complete isotropy, in the displacement correlation functions, 
is recovered at large scales. The translational correlation function is 



C 



Ko[r 



Ar 



(4.23) 



For the vortex lattice ( [2.4[ ), shear deformations dominate (cge -C Cn) in most of the phase 
diagram. The expressions for the function B^'^{r) which describe the crossover between the 
random manifold (intermediate distance) regime and the large distance regime then simplify: 



B 



L,Ti 



a ~. 
~ —h 



T,Li 



(4.24) 



In the limit of weak disorder ^ 3> a we find that there should be a well defined crossover 
function, i.e all curves should scale when plotted in units of x/^. The relative displacement 
correlation functions B^'^, as predicted by the variational method, are plotted in Fig. |T] and 
Fig. for the triangular lattice, by numerically integrating ( [4.19|) . 
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FIG. 1. Plot of for the triangular Abrikosov lattice, in the limit cge *C cn. The 

longitudinal (solid line) and transverse (dashed line) relative displacement correlation functions 
Bi t are defined in ( [4.91) . Kq is one of the first reciprocal lattice vectors, and ^ is the crossover 
length defined in ( [4.12| ). When x < one sees a power law with 2z^ = 1/3 characteristic of the 
random manifold regime. 
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FIG. 2. Plot of -y-Bl^t for the triangular Abrikosov lattice, in the limit cge <C cn. The 
longitudinal (solid line) and transverse (dashed line) relative displacement correlation functions 
Bl^t are defined in ( [4. 91) . Kq is one of the first reciprocal lattice vectors, and ^ is the crossover 
length defined in ( [4.12| ). When x > S, one sees the logarithmic regime. 



The crossover between the random manifold regime and the asymptotic quasi-ordered 
asymptotic regime is apparent, and occurs at a scale of order ^. At the length scale r = ^ 
where the random manifold regime cease to be valid, the translationnal order correlation 
function Cxir) = e~2^°'^p^°',i3(^\ is of order Ckq ~ 0.1. Therefore the crossover should 
be experimentally observable. In Fig. ^, we have shown the ratio R of the transverse to 
longitudinal displacements. 
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FIG. 3. Plot of the ratio Bt/Bl for the triangular Abrikosov lattice, in the limit cqq ^ cn, 
where Bl^t are defined in (|4.9| ). ^ is the crossover length defined in (|4.12| ). When x < ^ the ratio 
takes the random manifold value 2iy + 1. It decrease slowly to 1 in the asymptotic regime x ^ S,. 

As was shown by BMY@ii, its value is 2z/ + 1 in the random manifold regime (the vari- 
ational method give R = 4/3). At large scale, we find that this ratio decreases to i? = 1, 
and in that sense isotropy is restored. However, if one looks at the correlation functions 
for translational order, one finds that the difference between the longitudinal and trans- 
verse parts persists at large scales. Defining the longitudinal and transverse translational 
correlation functions by 

C^'^(r) = e-^^''"« (4.25) 

C^'^(r) correspond to correlation functions with a separation r, parallel and perpendicular 
respectively to the vector Kq. As is seen on Fig. ^, the ratio Rc = C^/C^ increases from 
one at short distances and saturates at a finite value at large distances. 
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FIG. 4. Plot of the ratio of the translationnal order correlation functions 

CL{r,z = 0)/C'r(r, 2; = 0) for the triangular Abrikosov lattice, in the limit cqq <C ch. Cl,t 
are defined in ( [4.25|) . ^ is the crossover length defined in (|4.12| ). The ratio increases from one at 
short distance x < ^ and saturates rapidly to a universal number in the asymptotic regime x 3> ^. 
The variational method predicts this number to be e. 



This value depends on the elastic constants, as seen from ( |4.23| ). In the limit cqq -C ch, 
this number takes a value which the variational method gives as athe universal constant e = 
2.7182... The fact that Rc saturates at large distance is a consequence of the existence of the 
quasi-order. Had a random manifold or a Larkin regime been valid up to large distances, this 
ratio would increase indefinitely. On the other hand if the system was genuinely ordered Rc 
would saturate at a much smaller value than e, a value which would go to one when T — > 0. 



As is discussed in more details in section |VD|, this should have observable experimental 
consequences. 



C. Physical Discussion: crossover lengths, dislocations 



One can give a simple physical explanation for the three regimes found here. Consider 
two flux lines separated in the ideal lattice by x. In presence of disorder the mean squared 
relative displacement is B{x). There is a length at which the potential seen by a line is 
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smooth. This length is the greatest of the correlation lengths of the random potential ^o, 
or the Lindemann length It = \J {u^)t- This defines a separation between vortices which we 

have denoted by / in this paper such that B{1) ~ max(^o^) ^t^)- At zero temperature it equal 
the length defined by Rc defined by Larkin Ovchinikov length ( and L^. in the z direction) and 
in general / can be thought of as the Larkin Ovchinikov length renormalized by temperature. 
Below this length the elastic manifold sees a smooth potential with well defined derivatives, 
thus a local random force can be defined. Indeed expanding in u the disorder potential energy 
in ( |2.18| ) gives a random force term f.u with f{x) = Y^x^{^)Kexp{—iKx) = W{Ri). In 
the sum over harmonics the maximum K is K^ax = 27r/^o- Thus this expansion is valid 
only as long as m <C .^o- This defines the range of validity of the Larkin regime, i.e at T = 
X < Rc and more generally x < I {we assume Rc> a). 

For separations larger than Rc but such that B{x) ^ each flux line explores only its 
immediate vicinity and feels different disorder. This is the regime explored by BMY which 
is identical to the random manifold. This can be seen, on a more mathematical level, from 
our model by summing over all the harmonics for instance on the replicated Hamiltonian 
( p.l3| ). One gets V{u) ~ J2r, ^iua ~ Ub — Ri). For u <^ a only the i? = term contributes 
and each line sees an independent random potential. This intermediate random manifold 
regime holds up to the length x = such that -8(0 ~ af which periodicity becomes 
important. There is no gain in energy to shift the lattice by a. In this regime displacements 
grow much more slowly and only the lowest harmonics contribute. This is the quasi ordered 
regime. 

In order to apply this theory to experimental systems one has in principle to worry 
about topological defects, such as dislocations. Although the influence of dislocations is 
still a controversial question, their influence has been clearly overestimated in the past. 
Let us mention some arguments, which we believe are incorrect, put forward to argue that 
unbound dislocations will proliferate even at weak disorder. An Imry Ma type argument is 
the following. The core energy cost of a dislocation cannot be avoided and scales as L'^~'^. 
A dislocation loop of size L creates extra-displacements of order 0(1) up to logarithms, in 
a region of size L'^. By adjusting the position of the loop one can hope to gain an energy 
from disorder L*^/^. Thus below d = A dislocation will be favorable. Such an argument is 
incorrect because it is again based on Larkin random force model for which disorder energy 
is linear in the displacement. For the real model the energy varies as cos(27rM/a) and adding 
a dislocation displacement will not necessarily gain enough disorder energy. In fact if the 
Larkin or the random manifold regime were true up to infinite scales it would indeed be 
favorable to create dislocations. The energy fluctuations due to disorder is ~ 2, '='-2+21/ ^ 
L'^~'^ . If z/ > dislocations will occur because it will always be energetically favorable to 
replace an elastic distortion by a dislocationil. However in the case of a lattice or if quasi 
long range order is preserved in the system as it is the case here, both energies scale the 
same way, since = 0, maybe up to logarithms. The prefactor of the disorder term can 
then be made arbitrarily small at weak disorder while the core energy of the dislocation is 
a given finite number. Thus if disorder is weak enough it is likely that dislocations will not 
appear. Even if they do the scale will be huge, and the effects associated with disorder that 
we discuss in this paper should be observable over a wide range of distances. If disorder is 
gaussian one could also argue that rare fluctuations will eventually lead to dislocations at 
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exponentially large scale. Realistic disorder however is bounded and thus such effect should 
be absent. 



D. Experimental consequences 

Let us now discuss in more details the experimental consequences of our findings. Two 
main types of experiments exist at the moment to probe the translational order of the 
vortex lattice: magnetic decoration experiments and neutron diffraction experiments. One 
would expect for these two experiments that the results of sections [IV B| and [IV C| apply. 



However, the direct comparison with experiments could be complicated due to the effects 
of the non-local elasticity and 3D anisotropy. These effects can be included, in principle, 
in the variational calculations by simply changing the elastic Hamiltonian, at the price of 
extremely tedious calculations. Even if a detailed treatment of such effects is beyond the 
scope of this paper, their importance can be estimated by the following simple arguments. 

In the high-Tc Abrikosov lattice, the elastic constants vary by orders of magnitude when 
the wave- vector goesfrom l/^o to l/-^- A good approximation of the elastic modulii for 

1 

47r 1 + A^g^ + K <1± 
cn(,) ^ (4,26) 



where A is the London penetration depth in the ab plane, Ac = LA along the c axis with 



r = ^JMz/M and c'^^iqz) is the single flux line contribution to the tilt modulus. One 
must have also B < H^ec = T'^^o/d'^ to avoid further effects of decoupling between planes, 
where d is the distance between CuO planes. Since in this regime Cge is dispersionless 
and much smaller than cn, most of the effects of non local elasticity comes from the 
momentum dependence of C44. In the region of Fourier space where C44 varies strongly, i.e 
for 1 -C Ac^x -C \c<i^ a good approximation to C44 is 

C44 ^ . fa 2 (4-27) 

this new momentum dependence of C44 will lead to very different lattice displacements. Using 
( p.l4[ ) with the proper C44(g_L) ( [4.27[) , one finds now [<j]{v) ~ f/log^(f). This leads, using 



( p.l7[ ), to a very slow growth of the relative displacements, {{[u{x) — ■u(O)]))^ oc log(x) or 
constant, for A < x < Ac. If the translational correlation length ^ > Ac, the random manifold 
regime will survive, whereas if ^ < Ac one would expect the non local elasticity effects to 
dominate the random manifold regime. Note however that the asymptotic large distance 
regime, for ^ > Ac will be completely unchanged. 

In decoration experiments however, one is usually in a regime of very small fields. In 
particular B < Hd (bulk), giving a > A. For example, in the Bell experiments^ performed 
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on lOyum thick samples of BSCCO, one has A = 0.3/im and Ac ~ 60A. The highest field 
picture (69 Gauss) have very large regions free of dislocations, for which one can hope to 
apply the theory of the present paper. For such fields a ~ 2A, and thus A < a < Ac- It is 
likely, since the interactions between vortices are less important than in the regime B > Hd, 
that single vortex contributions will dominate. Thus the q± dependence of C44 will be weaker, 
and one can hope non local elasticity to be unimportant. 

Recently, we have carefully reanalysedS the data of Ref.0, performing a Delaunay tri- 
angulation of the larger field images which do not contain dislocations. This allows to 
compute B{r) directly. Preliminary results indicate a very good fit to a power law be- 
haviour B{r) ~ r^*^ from r = a up to r = 30a with an exponent 2z/ = 0.4 ± .05. Thus, 
assuming dispersionless elastic constants, this is a strong indication that one is seeing the 
random manifold regime. The exact exponent 2i> for the random manifold regime is un- 
known, but the (Flory) value 2u = 0.33 predicted by the variational method is expected@ 
to be a (relatively good) lower bound. Another predictioniS for the exponent using refined 
scaling arguments, which might turn out to be more accurate, is 2i> = 4(4 — d)/{8 + m), i.e 
2z/ = 0.40 for (i = 3 and m = 2. The data exclude a Larkin type behaviour B{x) ~ x, and in 
fact there seems to be no measurable Larkin regime for small r, indicating that / ~ a. One 
does clearly observe a saturation in B{r) around r = 30 — 40a at a value of B{r) consistent 
with the predicted saturation to the slower logarithmic growth. However larger pictures 
would be necessary to conclude unambiguously on the crossover itself, as well as the large 
distance regime. The main obstacle is of a statistical nature, i.e there are not enough pairs 
of points uncorrelated statistically, to perform the necessary translational average. Larger 
pictures would allow such an average to be taken. 

Clearly both the understanding of the short distance regime, and of the importance or 
not of non local elasticity, and the existence of the quasi ordered regime deserve further 
studies. Other difficulties in interpreting data from decorations experiments can come from 
the fact that surface interactions may be different from the bulk onesS. It has been argued 
recently however that the effects of the surface interactions may be visible only at scales 
much larger than the size of the decoration pictures^. 

Another good probe of the correlations in the vortex lattice, which is free of potential 
surface problems, is the neutron scattering experiments. Detailed neutron diffraction studies 
are now available for different type II superconductors, such as NbSe2 (see reference |3^) as 
well as BSCCO. Neutron experiments measure (up to a form factor taking into account the 
field distribution created by a single vortex line) the Fourier transform aX k = Kq + g of 
the density correlation around a reciprocal lattice vector Kq. The structure factor which is 
measured is given by 

S{q) = J rf3^gigxg-ix„x^B,^(x-) (^4 28) 



where B is given in (|4.9| [4.18| |4.19| |4.20|) . The full calculation of S{q) requires a numerical 



integration of ( [4.28 ), but the main features can be given analytically. Let us recall that 
^ is the translational correlation length due to disorder defined in ( [4.12|) and that trivial 
anisotropy has been taken into account by proper rescalings of z versus r directions. 

At small q, q < 1/^ the integral ( |4.28| ) is dominated by the large distance regime where 
B{x) = A^logi^x), where A3 = 1 according to the variational method. The structure factor 
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is therefore 



Siq) ~ (1/g) 



3-A3 



(4.29) 



and thus diverges at small g, a consequence of the persistence of quasi-long range order 
in the system. True Bragg peaks therefore exist. This is in sharp contrast with previous 
predictions assuming simple or stretched-exponential decay of the translational correlation 
function up to large distance. At a wavector of order q ~ 1/^, the behavior of S{q) will 
crossover to a slower decay, controlled by the random manifold regime. In this regime 



A clear signature that one is indeed in the regime ( [4.29| ) described here should show in the 
neutron experiments by the fact that S{q) has no true half-width. The maximum value 
of S{q) will be limited either by the experimental resolution of by the distance between 
unpaired dislocations ^n. Varying these parameters should leave the rest of the curve nearly 
unchanged (this is valid as long a.s C,£) C,). The distance between dislocations could be 
controlled by annealing the lattice, either using a driving force or a field-cooling procedure 
similarly to what is done in Ref. |3^. 

Another interesting prediction can be made for the in-plane, g_L-dependent ratio S{q± \\ 
KQ)/S{q± -L Kq). After integration over q^ the structure factor becomes, at small q± 



S{q) ~ (1/g) 



3+2i/ 



(4.30) 




(4.31) 



Thus the ratio goes for small q± to the value 



S{q^±Ko,z = 0) 
S{q^\\Ko,z = 0) 



= e 



(4.32) 
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V. FUNCTIONAL RENORMALIZATION GROUP 



Another method widely used to study disordered problems is the functional renormal- 
ization group. It turns out that its application to the present problem due to the periodicity 
in (|2.13| ) is simple. It provided a good complement to the variational method, none of the 
methods being rigorous. The functional RG can only give results in an e = 4 — ci expansions, 
which does not have presently the status of rigor of the standard e expansions of usual critical 
phenomena for pure systems. In particular the effects of multiple minima will affect higher 
orders in perturbation theory and could very well result in replica symmetry breaking insta- 
bility in the FRG flow, as found recently in Ref.0 (see also Section |V|). On the other hand 
the functional RG should include fluctuations more accurately than the variational method, 
provided it does not miss another part of the physics. Comparison of the two methods near 
four dimensions should allow to test their accuracy. 



A. One component model 

For simplicity we confine our study to a model with isotropic elasticity as in (|2.13| ). Let 



us first consider u to be a scalar field (m = 1) and set c = 1. The full replicated Hamiltonian 
is 

H/T = ^J rf^x(V«(x))2 'W^^l ^'^^("'^(^) - "^(^)) (5-1) 

ab 

For simplicity we take Kq = 27t, so that the function A{z) is periodic of period 1. In the 
original Hamiltonian A{z) is a sum of cosine given in (|2.13| ). One then performs the standard 



rescaling x — e'x and u e^'w. The idea of this renormalization is to perform an expansion 
around a classical solution at zero temperature. One should keep the whole function A in 
the renormalization procedure. The RG equations to order e = 4 — d have been derived by 
D.S. Fisheiii for the random manifold problem (see alsoc3cZI)Ea 

^ = (e - 4C)A + C^A' + l-{A"y - A"A"(0) (5.2) 

A factor 1/5^ = 1'^-^'K^I'^V\dl1\ has been absorbed into A in (|5.2| ). The temperature is an 
irrelevant variable and flows to zero. The correlation function 

f (g) = TG(g) = {u\q)u\-q)) (5.3) 

satisfies the RG flow equation 

f(g, T, A) = e(^+2C)if (ggi^ Te^""'^)', A(/)) (5.4) 

The periodicity of the function A implies that the roughening exponent is C = for 
the large distance behavior. This allows us to obtain the only periodic fixed point function 
A*(^) in the interval [0, 1]: 
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^*{z) = ^{^^-z\l-zf) (5.5) 

Values for other z are obtained by periodicity. The fixed point is stable except for a con- 
stant shift, which corresponds to a change in the free energy. The linearized spectrum is 
discrete and the eigenvectors can be obtained using hypergeometric functions. The fixed 
point function is non analytic at the origin. It has a singular part which behaves as z^\z\, 
for small z and leads to A**-^''(0) = oo. As discussed below this is a general feature of fixed 
points for this type of disordered systems^. For a periodic fixed point, i.e. C = 0? oiis can 
set e'*g = 1/a in ( p.4| ) this allows to obtain perturbatively, provided I* ^ 1 

f(g,T,A)= f(i,T^O,A*) (5.6) 

y gay a 



One can evaluate the correlation in ( |5.6|) , at a scale of the order of the cutoff, perturbatively 
in A. One then gets 



r(g,r,A)= - 



1 Y -A"*(0) 



qa I 



(5.7) 



k=l/a 



Using (|5.5| ), and remembering the factor 1/5*^ in A, one obtains 



This gives 



((M(x)-n(0))2) = -log|a;| (5.9) 
whereas the variational method gives 



{iu{x)-u{OW) = —\og\x\ (5.10) 



( ^.9|) gives a power law decay for the translational correlation functions with an exponent 



Ad,RG = e(27r)^/36 = l.lOe against A^/jg = ^ for the variational method. The agreement 
of the two methods on the exponent is within 10%, which is very satisfactory. The fact 
that AdyAR < Ad^RG is not surprising since the variational method underestimates a priori 
the effect of fiuctuations. One can remark that omitting the term (A")^/2 in ( p.2|) leads 
to a fixed point A* (2;) = cos(27rz)/(27r)^, which gives exactly the same exponent than the 
variational method. 

At intermediate distance it is enough to focus on the small u behavior of the function 
A, and thus to forget in effect the periodicity. At short distances the function A is analytic. 
In that case@ 

^^ = (e-2C)A"(0) (5.11) 
Setting C = e/2 allows to get a fixed point A (2;) = Az"^ — 2A^. Using ( p.4[ ) one obtains 
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f (,) ~ ^ (5.12, 

which corresponds to the Larkin random force regime. This however holds only at short 
scales. This fixed point is unstable, and a non analyticity at z = develops, corresponding 
to an algebraic decay of the Ak in (|2.13| ). A eventually flows towards the long distance 
regime described by the flxed point (|5.5| ). There might be an intermediate random manifold 
regime. 

Another renormalization method that has been used, was a real space RG by Villain and 
Fernandezi^. For 2 < c? < 4 this method, which is approximate, also predicts a logarithmic 
growth of the correlations. It does not allow however to compute the universal prefactor 
Ad or the crossover function. The agreement between these methods, none being rigorous, 
lends credibility to the additional results in d = 3 obtained using the variational method. 



B. General case 



Let us consider now the more general case of an m component vector u, and isotropic 
elasticity (|2.2|) . The equation giving the flxed point function becomes instead of ( |5.2| ) 



A[u] + ^{d^dpA[u]y - d^df,A[u = OWpA[u] = (5.13) 



while the displacement correlation function becomes,similarly to ( ^.7|) 



{u^{q)up{-q)) = {d^df,A*[u = 0])^ (5.14) 



where e has been included in A. For the case m = 1, (|5.13| ) reduces to ( ^.2| ). For the case 



m = 2, the analysis depends on the symmetry of the lattice. For a square lattice a separable 
function 

A[u^,Uy] = A*[u^]+A*[uy] (5.15) 



where A* is a solution of ( ^.2|) satisfles ( ^.13|) and gives (|3.19| ) with the same exponent 



Ad^RG than for the case m = 1. A rectangular lattice would give the same exponent. The 
triangular lattice is more difficult to treat and no simple solution of ( p.l3| ) can be found. We 



have performed a numerical solution of (^.13 ). There is a non trivial solution which has the 



full symmetry of the triangular lattice. In that case one would expect in general different 
exponent A^ than for the square lattice, unless there is an hidden symmetry reason, for 
which the exponent does not depend on the lattice symmetry. It is difficult to get a high 
precision for the exponent because of the non analytic nature of the solution. The numerical 
value found for Ad was within 5% of the one for the square lattice, but we were unable to 
decide within our accuracy whether the two exponents were equal or different. The exponent 
is again very close from the one predicted by the variational method, which is independent of 
the lattice symmetry. Once again, neglecting the term ^{dadj3A[u])'^ in ( p.l3| ), one recovers 
exactly the result of the variational method. 
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VI. D=2 



In (i = 2, thermal fluctuations are expected to play a more important role. Already in the 
case of the pure system, they change the true long-range order of the lattice into a power 
law decay of the correlation functions, with an exponent controlled by the temperature. 
One can therefore expect a stronger competition between disorder and temperature than in 
higher dimension. In addition standard renormalization group techniques are available in 
d = 2 and can be compared with the variational method. In the section |V1 A| we examine 
the d = 2 problem using both the variational method and the renormalization group. We 
will focus mainly on c? = 1 + 1 (flux lines in a plane). The results are mostly relevant there 
since the starting model ( |2.13|) becomes exact, due to the fact that dislocations cannot exist 
in d = 1 + 1. The physical consequences for various experimental systems both in ci = 1 + 1 
and d = 2 + will be discussed in section |VI B|, together with the effects of dislocations. 



A. Theoretical results 

In d = 2 the variational method applied to the starting model ( |2.13 | 1 leads to a solution 
which belongs to the class of "one-step" replica symmetry breakingSS, in some extended 
sense, i.e such that [cr]{v) = vanishes for v < Vc and [cr]{v) > for t>c < f < 1. This 
can be seen readily by taking the limit (i — 2"*" in (|3.15|) , a limit which vanishes for v < 
vo = TKq/{4:7ic). This solution represents a glass phase. Since vo cannot be larger than 
1, the glass phase exists in d = 2 for T < Tc = Attc/Kq. For T > Tc the disorder is 
irrelevant, and the replica symmetric solution is stable, as already discussed in section [Til B 



The detailed behavior below will again depend on whether one considers simply a single 
cosine model, or takes into account all the harmonics present in ( p.l3|) . For simplicity we 
will focus here on the single cosine model, which has been simulated numerically and is 
interesting in itself. The main effect of the higher harmonics is again to allow for a random 
manifold crossover regime, at low enough temperature in the glass phase. It is examined in 
details in appendix 0. 

In the case of the single cosine model ( ^.llj ) one can look simply for a constant solution 



[(T](f ) = Si for V > Vc- The details of the calculations can be found in appendix |C|, where 
the saddle point equations ( PU[ ) for Si and Vc are derived. These equations are solved in 
( |CT2|) to give 

^ ° mT vsi y 

where the cutoff A ~ c(27r/a)^ Assuming the cutoff to be very large (i.e ^ <^ a) in ( |6.1| ) one 
gets 

- - ° (6.2) 



47rc 

1/(1-T/T,) 

/ / \ /A ^ \ 

Si = A 



I 47rmA ] 
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As we will show below Si defines a characteristic length 



e(T) = = a(c2aVA)i/(2-TX„V2.c) ^ a{^a^ / Af'^^-^'^'^^^ (6.3) 

above which there is a crossover to the asymptotic regime dominated by the disorder. The 
above expression for ^(T) in o? = 2 is valid for ^ ^ a. At zero temperature it coincides 
with the length found using simple Fukuyama Lee arguments (see section [II C|) and it is 
renormalized downwards by thermal fluctuations at finite temperature, an effect specific to 
two dimensions. 

Using (|2.14|) and the form (|C2D for G{q) one obtains for the relative displacement 
1 



Bix) = —{(uix) -uiQW) (6.4) 
m 

where Bq{x) is the value of B in the absence of disorder 

5o(a;) = — log(Ax) (6.5) 
cvr 

( |6.4|) is convergent although each term is individually divergent, but is easily regularized by 
multiplying by Jo{ek) with e — ^ 0. This leads to 

B{x) = — \og{Ax) + (log(a;/0 + Mx/O + 7 + log(l/2)) (6.6) 

cvr cvr 

where 7 is the Euler constant. The expression ( |6.6|) gives the following crossover for B{x) 

T 

B{x) ~ — log(x) x < ^ (6.7) 
cvr 

Tc 

B{x) ~ — log(x) X 3> ^ 
cvr 

where Tc = Attc/K^. Note that for a ^ x ^ ^ there is in principle a Larkin regime where one 
has algebraic growth of the disorder part of the correlation function with 2i>l = {4 — d) = 2 
plus logarithmic corrections: 

Bix) = B,{x) + ^^(log(e/(2x)) + 1 - 7)(|)' (6.8) 

However, except at very low temperatures, the thermal part always exceed the disorder part 
and thus disorder effects are masked by thermal effects at short distances. 

The variational method predicts therefore a simple logarithmic growth of the displace- 
ments at large distances, both above Tc and below. The effects of disorder are limited to 
the freezing of the prefactor to the value Tc for temperatures below Tc- Note that Tc is a 
universal quantity, independent of the strength of the disorder. Such a result is valid in the 
limit where the ultra- voilet cutoff A is very large. The disorder strength enters the crossover 
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length ^ above which the asymptotic behavior for T <Tc can be observed. Of course ^ — oo 



when T — >■ Tc as can be seen from (6.3). 

Note that the effect of the cutoff, which could be important for a numerical simulation not 
at small disorder, lead to some temperature dependence of the amplitude of the logarithm. 
The amplitude in ( |6.7D becomes 

A = T/iv^cTT) = TJicTT) .fy^,,. (6.9) 

log(l + Si/A) 

where one can use Si from ( |6.2| ) or better from the equation ( |(J13| ) for y in the appendix y. 
One finds an increase of the amplitude A when the temperature decreases. 

In c? = 1 + 1, it is also possible to write renormalization group equations for the 
disorderi§'00. To lowest order, and on a square lattice for simplicity, such RG equations 
were derived by Cardy and OstlundS and read: 

§.,2_|I)A_CA^ (6,0) 

where we denote A = Akq, and 



dl T2 



(6.11) 



Both A and Aq are defined in (|2.13|) . C is a constants which is unimportant for our 
purposes. For T > the disorder is irrelevant, in agreement with the variational method. 
For T < Tc the disorder term becomes relevant and there is a new non-trivial fixed point 
at A = (Tc — T). This fixed point however has the unusual feature that the variable Aq 
flows to infinity Ao(/) oc /. However since this variable does not feed back at any order in 
perturbation theory (only averages of the type (J2a Ca(t>aY with C'a = appear) it has 
been assumed that this fixed point was correct. 

Using the RG one can again define a short and large distance regime. At short distances 
the RG is certainly correct, and is more accurate than the variational method since it 
treats the fluctuations correctly. As was noted in section [111 C 1| , at short scales x such 
that \u{x) — m(0)| -C ^ (for the single cosine model), it is possible to expand the cosine 
to equivalently recover the Larkin random force model. The correlation function for that 
model reads simply 

G(g)(x4 + 4 (6-12) 

This usually leads to SB{x) ~ x'^, but here one must take into account renormalization by 
thermal fluctuations. This is done by integrating the RG equation in the small distance 
regime where we note that (|6.10|) is correct as long as Aa^/T^ <^ 1, which is equivalent to 



{a/^Y /(iT/aY -C 1. One can integrate ( |6.10| ) to obtain 



A(/) = A(0)e'(2~2^) (6.13) 

Applying the RG flow equation 
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G{q, T, Ao, A) = e2^G(ge^ T, Ao{l), A{1)) (6.14) 
where G has been defined in ( p.l4|) , immediately leads to 

SG(,) ~ Ml ^ _i_ (6.15) 

The RG therefore predicts that the Larkin regime is in fact anomalous with an exponent 
continuously varying as a function of the temperature 

513{x) = B{x) - Bo{x) ~ a2(x/0^"^ (6.16) 

instead of (|6.8| ), for x < ^ (in the low-T regime ^ is replaced by another length). There are 
also corrections coming from the renormalization of Aq. Integrating (|6.11|) one gets 



i^2a^A(0)_2 

T ■ 



Ao(0 = Ao(0) + ^y^'^^-'l) ^"''^'"^^ ~ ^^-^^^ 



but such corrections are obviously smaller at short distances. 

At large distance Ao(/) 3> 1, in order to obtain the correlation functions, one has to 
assume that the unusual CO fixed point is indeed correct. If one does so, correlation functions 
be compute d§ using RG fiow equation (|6.14| ). Iterating until /* such that e'*g = 1/a 



can 

allows to obtain for large /: 



G{q, T, Ao, A) = \G{a, T, Ao(/), A*) ~ log{l/q)/q' (6.18) 



which leads to B(x) ~ log (x). In ( p.l8| ) it has been assumed that simple perturbation 



theory could be done for the correlation functions at scale x = a. The RG approach would 
therefore predict a log^(x) growth of the displacements, at variance with the predictions 
of the variational method which gives a simple log. In fact the RG results is based on 
the assumption of replica symmetry. As we have shown recently^, a careful analysis of 
the Cardy Ostlund fixed point and of the RG fiow shows that it is unstable to replica 
symmetry breaking. When RSB is allowed one obtains a runaway fiow of the RG which is 
consistent with the findings of the variational method. Two recent numerical calculations 
on this model0'ii seem to confirm that the GVM does describe the correct physics at large 



distance. None of them is compatible with a log (x) growth of the displacements. In Ref. |5T 
no change in the static correlation functions was observed in the presence of small disorder, 
whereas a transition occurring in the dynamic correlation functions was observed at Tc. 
A careful comparison was performed with the predictions of the RG calculations, and the 
results were found incompatible. These numerical results are however consistent with the 
prediction of the variational method. Indeed for such weak disorder the length ^(T) is very 
large, especially near Tc, and simulations performed on a too small system will show no 
deviations as is obvious from (|6.7| ). However in Ref. ^ the disorder is much larger, and 
^(T = 0) is of the order of the lattice spacing a. This simulation indeed shows quite clearly 
a freeze of the amplitude of the logarithm below at the value A = Tc/{ctt). 
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B. Physical realization in 2 + dimensions: Magnetic bubbles 



Elastic models in two dimensions have been studied for some time. The Hamiltonian 
( p.l3| ) describes several physical disordered systems, such as randomly pinned flux arrays 
in a planeiiiii, the surface of crystals with quenched bulk or substrate disorder^, planar 
Josephson junctional, domain walls in incommensurate solids. In addition, a very nice re- 
alization, investigated in detail recently, is provided by magnetic bubbleJil. In such systems 
one should be able to test the predictions of the previous section. However, if the elastic 
objects are not lines {d = 1 + 1) but points (d = 2 + 0) it is now important a priori to take 
into account the effects of free dislocations. In two dimensions, dislocations are expected 
to be much more important than in three dimensional systems, and to have observable 
consequences on the destruction of translational and/or orientational order. A very naive 
argument in d = 2 is that while the energy of a dislocation pair of separation r increases as 
~ TQLog{r/a) in the absence of disorder (a slowly growing function but still allowing for a 
quasi-solid phase for T < Tq), in the presence of disorder, the energy should saturate to a 
finite value -Emax — ToLog^^/a) when r > This is because long range order is supposed to 
be destroyed beyond this length, effectively screening the elastic pair interaction. If it is the 
case, unpaired dislocations will be thermally excited at any finite temperature with weight 
Q-Emax/T^ Note that the same type of argument in d = 3 excludes dislocation loops of size 
r whose energy scale as rLog{^/a) and thus cannot proliferate so easily. This was discussed 
in section liV C. 



As was discussed in section |VI A| , in the absence of dislocations quasi-long range order 
persists in the system even in the presence of disorder, and the above naive arguments need to 
be reexamined. A quantitative theory including both disorder and dislocations is difficult. 
One can however get an idea of the effects of dislocations by using the renormalization 
equations of Cardy and Ostlund in the presence of topological defectst^. These equations 
were derived by Cardy and Ostlund for a one component field (XY model) with both disorder 
(random field) and vortices. They correspond to the modification of ( |6.1CI| ) and ( |6.11| ) to 
include the dislocation fugacity, and an additional equation for the renormalization of the 
dislocation fugacity itsellp3. Generalizing them to a triangular lattice beyond linear order 
goes beyond the scope of this paper, but we do not expect radically new conclusions, rather 
a change by a small factor of the parameters ( temperature, etc..) (note that extension 
to a n = 2 components on a square lattice is obvious). There are several reasons why we 
believe that the CO equations might overestimate the effect of dislocations. One is that, 
as we saw before ( |VI A| ) these RG equations implicitely assume replica symmetry. They 
lead to a faster decay of the translational correlation function in the absence of dislocations, 
of the type exp(— log^(r)) and thus overestimate the effect of disorder compared to the 
GVM, which includes replica symmetry breaking. The GVM predicts that quasi long range 
order exists in the system. Thus the CO equations should also overestimate the effects of 
dislocations. 

A stability diagram can be constructed by examining the renormalization equations to 
linear order, for the fugacity of topological defects and for the disorder. In the system 
considered by Cardy and Ostlundil the fugacity of vortices y satisfies to lowest order in y: 
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-dl = - 2KlT ^ 2K!f^^ ^'-^'^ 

while the g = Kq component of the disorder Akq renormahzes to hnear order as: 

^^2A.„,1-MI, (,.0) 

The influence of internal disorder on a two dimensional crystal was studied in Ref.0. 
In that work the disorder was taken as a quenched random stress coupling to the strain, 
i.e a term —a{x)'V ■ u was added to the hamiltonian with cr(x)o"(x') = po'^Ao6{x — x'). As 
discussed in this corresponds to include only the long wavelength part of the disorder, 
and to neglect the q = Kq component of the disorder, which is very important for the 
present problem of substrate disorder. For the 2D elastic triangular lattice the fugacity of 
dislocations y renormalizes, to linear order, in a way very similar to |6.19 : 



-2y(l C66(cii-C66) ^ pgApg^ C66 2 

dl AnT cii ttT^ ch 

One easily sees that, to linear order, the q = Kq disorder renormalizes also very similarly to 

The stability diagram is shown in figure ^. 
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FIG. 5. Stability diagram of a two dimensional solid with weak quenched substrate disorder, 
as a function of the temperature T and the long wavelength (i.e harmonic) part of the disorder 
Aq. The diagram indicates schematically the stability of the harmonic part of the free energy to 
dislocations (fugacity y) and short- wavelength disorder (y ~ ^Ko) by showing the relevance of 
these variables to linear order. 



In presence of dislocations there are thus various possible regimes. There are two critical 
temperatures. One is the KTHNY melting temperature = o^cqq^cu — C66)/47rcii, above 
which dislocations unbind for the pure system which melt into a hexatic. The other one is 
the glass transition temperature Tq = ShcuCqq/ Kq{cu + cqq), below which disorder becomes 
relevant in the absence of dislocations. The ratio between these temperatures is equal to 
half of the exponent 77^0 of the pure system at T^. For the CO model it is thus universal and 
equal to T^/Tg = 1/8. For the triangular lattice it depends on the (renormalized) elastic 
constantsEj but cannot exceed 1/6. In figure |^ we have shown schematically how y and 
y ~ Akq are renormalized, at linear order as a function of temperature and long wavelength 
disorder Aq from 6.21 3.22 . At high temperature dislocations are relevant and disorder 
is washed away by thermal fluctuations. Below Tg disorder becomes relevant leading, in 
the absence of dislocations, to the Cardy-Ostlund line of fixed points. This line however 
is unstable to dislocations. In the shaded region in figure ^ the dislocations would be 
perturbatively irrelevant for the system with A^^^ = and the solid would survive ( in fact 
in that case if one includes non linearities in y the solid is stable in an even greater region@ 
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). However in this region the disorder y ~ Ak^^ is relevant. It will eventually increase Aq as 
can be seen from ( |6.11 ), and drive the flow towards the region where y itself increases and 



unpaired dislocations will eventually appear at large scale. 

Let us now estimate at which scale dislocations will appear, assuming weak disorder. 
The disorder = becomes of order one at lengthscale of order ^. The key point is that 
up to this lengthscale the fugacity of dislocation has been renormalized downwards and is 
now much smaller with y(^)/y{a) ~ {a/^Y^^~'^"^^'^\ At this lengthscale, one ends up with a 
system for which the disorder is of order one and the fugacity of dislocations is extremely 
small. One can therefore predict that the typical distance between unpaired dislocation 
in the shaded region is much larger than ^, the distance above which effects of disorder 
manifest and ( |6^ can be observed. It is impossible to compute rigorously the ratio ^d/C 
using perturbative RG since beyond ^ one of the coupling constant (Akq) is large, but one 
can still estimate this ratio by the following heuristic argument. If one assumes that, when 
has become of order one, ( |6.11| ) is still valid, then above length of order ^, 



Ao(/) = Ao + a/ (6.23) 

where a, is a coefficient of order one. Then estimating the distance between dislocations as 
the scale at which y becomes of order one, one gets using |6.19| 

~ eexp (C W/2((e/a)2(^-/^-^Vz/(«))) (6-24) 

where C is a constant of order unity and y{a) the bare value of the dislocation fugacity. 

The resulting prediction is that in the shaded region of figure ^ dislocations appear only 
at scales large compared to ^, and the main reason for decay of translational correlations is 
the pinning of the elastic manifold. Thus the elastic theory developped in this paper should 
apply up to distances up to order C,d- In that low temperature region perturbative RG does 
not allow to obtain precisely the correlation functions even at not too large distance since 
one is far from the perturbative region T ~ T^. The variational method on the other hand 
predicts the following results: as shown in appendix ^ there is a crossover temperature 
T* = Tg / log{C, / a) . For temperatures T > T* there is no intermediate random manifold 
regime and the displacement correlation function should show the logarithmic growth of 
(|6.7|). However since the disordered solid regime corresponds to T < — Tg/8, one is 
likely to be in the regime T < T*. In that case there is a random manifold regime at short 
distances where the displacement correlation function grows as B{x) ~ (x/^)^/^. At large 
distance one recovers the logarithmic asymptotic regime. The anisotropy ratio between the 
longitudinal and transverse displacements R = Bt/Bl should crossover from 2z/ + 1 ~ 1.7 
in the random manifold regime to 1 in the asymptotic regime. 

Experimentally, in the system of magnetic bubblesilH some regimes are found indeed 
where the distance between dislocations£i:) is of the order of up to five times the translational 
correlation length. The authors of Ref.llj observe a crossover between a regime where ~ C 
to a regime where C,d ^ C,, which probably corresponds to the transition to the shaded 
region discussed above. Indeed, the measured ratio Cd/C keeps increasing with density and 
is limited there only by experimental setting. It would be very interesting to increase the 



range experimentally accessible. In Ref. 14 some comparison with the random manifold 



result was also performed. However the definition of the anisotropy ratio used in Ref. P 
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(Ct/Cl) does not seem appropriate, since the data were first fitted to a simple Lorentzian 
shape. New insight in the physics of such disordered 2D systems could be gained by a 
reanalysis of these data. More theoretical work is also needed since the above analysis of 
the dislocations is very crude. It has been generally assumed that the non perturbative RG 
flow in presence of dislocations was towards a kind of hexatic, but nothing is really known 
on this phase. In particular, in the Cardy-Ostlund approach, the vortices are supposed to 
be completely thermalized and to see only a uniform potential, and feel no effect of the 
disorder. This is physically incorrect, since in a disordered system the fugacity will also 
depend on the position. One should therefore take also into account terms which will pin 
the vortices. This is another reason why we believe that the CO analysis overestimates the 
effect of topological defects. 

To conclude, d = 2 elastic disordered systems have therefore a very rich behavior. Many 
other experimental systems could be investigated, among them coUoidsBil and thin super- 
conductors. A realization which would exclude dislocations could be to adsorb a polymerized 
membrane, with very high dislocation core energy, on a disordered substrate. 
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VII. CONCLUSION 



In this paper we have developed a quantitative description of the static properties of a 
lattice in presence of weak disorder. We have derived a model (p.l3| ) valid in the elastic 
limit Vn ~ a/S, ^ 1 which contains the proper physics at all lengthscales, and allows 
to describe the various regimes as function of distance. We have apphed the Gaussian 
Variational Method (GVM) to this model and computed the correlation functions of the 
relative displacements. This method has the advantage of being applicable in any dimension. 
The comparison with the renormalization group study that we have also performed, which 
is possible in d = 4 — e and in c? = 2, indicates that the GVM is an accurate variational 
ansatz for this problem. It also shows that the GVM should be a good tool to explore 
other disordered problems where different length scales are present. Contrary to previous 
studies the present analysis includes both the effect of metastable states and of the intrinsic 
periodicity of the lattice, both effects being found to be very important. 

We found that the effect of impurities on the translational order of the lattice is weaker 
than was previously thought and that quasi-long range order persists at large scales. The 
resulting phase, which we call a "Bragg glass" has both the properties of being a glass and 
of being quasi-ordered. The analogy with the quasi-order which subsist at low temperature 
in pure two dimensional solids is puzzling. The Bragg glass possess several intrinsic length 
scales. At very short scales it behaves as predicted by Larkin's random force model. This 
regime should be quite limited in (i = 3 at low temperature and for a potential rough at scales 
smaller than the lattice. In c? = 2 however, thermal effects being stronger, this regime is 
wider. At intermediate scales, and low temperature in ci = 2, the system behaves as a random 
manifold of flux lines independently pinned by impurities with stretched exponential decay 
of translational order. At distances larger than ^, we have shown that due to the periodicity 
of the lattice, the pinning by impurities becomes less effective and quasi order survives. 

We have not attempted to include explicitly topological defects in the present quantita- 
tive study. However, we have given evidence both in = 3 and d = 2 that for weak impurity 
disorder the effects of dislocations is less important than is usually beheved. This is due to 
the fact that quasi-long range order subsists even in the presence of disorder for the elastic 
system. In c? = 2 it is possible, but by no means established, that unpaired dislocations 
will appear at large scale. We have estimated conservatively the length scale at which this 
happens and found that it can be much larger than ^. Thus in ci = 2 there is a regime 
where the decay of translational order is due primarily to the elastic displacements induced 
by disorder. In ci = 3, we have argued that there is probably a phase without unpaired 
dislocations for weak disorder and at low temperature. 

Thus in the three-dimensional high temperature superconductors there could be gener- 
ically be (at least) two types of glass phases caused by point impurities, i.e two types of 
vortex glasses. The first one is the strong disorder vortex glass phase which is described by 
the gauge glass type of mo delffl and presumably contains a lot of dislocations. The second 
one is the quasi-ordered Bragg glass described in the present paper. A natural speculation 
then is that the tricritical point recently observed in experiment^, whose position in the 
H — T plane can be raised by adding more point impurities in a controlled fashionlil, has 
something to do with the relevance of dislocations. It could mark the separation between 
these two phases since it is natural that the quasi-ordered glass melts through a sharper 
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transitioitJ, much like a pure solid, while the vortex glass transition of Ref.QuJ is believed 
to be continuousEH In Lt, it is not even clear whether, strictly speaking, the strongly 
disordered vortex glass phase of Ref.i'0 is a true thermodynamic phase, rather than a very 
long crossover from the flux liquid. Indeed recent simulations indicate that d = 3 could be 
slightly below its the lower critical dimension^. The Bragg glass, on the other hand, should 
not suffer from such existential problems. 

The results of the present study suggest some comments concerning experiments. First 
decoration experiments should be reanalyzed. Clearly a fit to a simple exponential as in- 
spired by Ref.S0 is certainly inadequate for pictures where no dislocation is present. Since 
we predict that the crossover to the asymptotic quasi-ordered regime occurs when the cor- 
relation function CKo{r) is of order 0.1 this crossover should be observable, in principle. 
One must keep in mind that Abrikosov lattices might not be the ideal system to compare 
with the theory, since they are very complicated objects with additional intrinsic scales. In 
particular the effect of nonlocal elasticity, 2d to 3d crossover effects and surface interactions, 
although they can in principle be incorporated into the method, complicate the analysis by 
introducing new crossover lengths such as the London length. It would be therefore highly 
desirable, in order to test the above predictions to investigate simpler elastic systems. One 
example are colloids which have been very successful to investigate translational order in 
pure systemsim. In addition, these experiments could help decide whether dislocations 
are a thermodynamic property of the state or simply a non equilibrium feature which can 
be eliminated. Neutron diffraction experiments could also help decide on these issues. We 
obtained several characteristic predictions for neutron diffraction, such as the existence of 
algebraically diverging Bragg peaks and predictions for the ratio of transverse to longitudinal 
scattering intensities. 

We have addressed here only the statics of a lattice in presence of disorder, the quantities 
of interest in that case being mainly related to positional order. It would be clearly very 
interesting to obtain information on the dynamics and especially the driven dynamics of such 
systems. This is indeed particularly important for superconductors where most experiments 
measure only dynamical quantities. Several remarks are in order. 

There exist presently a qualitative approach to the dynamics of the vortex glass based 
mostly on scaling arguments for energy barrierdl. On such a qualitative level one usually 
assumes a single scale for the energy landscape, i.e that barriers Eb{L) ~ scale with the 
same exponent 6 which describes energy fluctuations in the statics. Assuming that there is a 
regime of transport where plastic deformations can be neglected, i.e that energy barriers are 
controlled only by elastic motion, the present study of the statics can be used to describe the 
creep regime at low temperature. The argument goes that in presence of an applied external 
force /, such as the Lorentz force created by an external current |/| ~ \j\, the barrier for an 
optimal deformation u ^ L'^ away from a low energy configuration on a scale L is lowered 
and becomes Ei,{L) ~ fuL'^. To unpin locally the manifold thus needs to go over an 

Arrhenius barrier Eiy ~ 1//^ with = 6/{d + u — 6). This, and the exponents found here 
lead to a non linear voltage-current relation V ~ exp[— l/(Tj^)] in the flux creep regime, 
where fi crosses over from fi ^ 0.7 — 0.8 to /i = 1/2 as j decreases. 

On a more fundamental level one anticipates that the methods used here for the statics 
and the results thev yielded will set the ground for the study of the dynamics. As was 
discovered recently!^ the hierarchical structure of states which is encoded in the replica 
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symmetry breaking solution of the statics finds an exact translation in the dynamics. There 
it corresponds to the existence of a long time dynamics for which the fluctuation dissipation 
theorem breaks down and non equal time correlation and response functions become highly 
non trivial. It would be interesting to generalize these studies to the flux lattice. Let us finally 
note that the quantities computed here using the statics are equal time disorder-averages and 
they will coincide with the observed translational averages in a given experiment provided 
equilibrium has been reached, which will always hold below a certain length scale. 

We thank D.J. Bishop, J. P. Bouchaud, M. Charalambous, D. Fisher, M. Gabay, D. Huse, 
T. Hwa, M. Mezard, CM. Murray, R. Sheshadri. for interesting discussions 

APPENDIX A: RELABELING OF THE LINES 

In this appendix we detail the derivation of expression (|2.^ ) for the decomposition of the 
density in term of the relabeling field ( |2.7| ). We denote d-dimensional positions x = {r,z) 
where r belongs to the m-dimensional transverse space. The density is given by (|2.3| ) 

p{x) = J26{r-R,-u,{z)) (Al) 

i 

In order to take the continuum limit, one can introduce a smooth displacement field u{r, z) 
by 

= L (S^^"" ^ e--«^n,(z) (A2) 

such that u{Ri,z) = Ui{z) which has no Fourier components outside of the Brillouin zone 
(BZ). In terms of the smooth field one can introduce the relabeling field 

(f){r, z) = r — u{(j){r, z), z) (A3) 

In the absence of dislocations there is a unique solution of ( |A3| ) giving u{r, z) as a function 
of (p{r^z). is a m-component smooth vector field labeling the lines, and which takes an 
integer-like values at each location of a line 

(t){Ri + u{Ri,z),z) = Ri (A4) 

Substituting ( |A^ ) in ( [ATI) one gets 

p(x) = - ^))det[9„0^(r, z)] (A5) 

i 

Using the integral representation of the 5 function, ( [A^ ) becomes 

p(x) = dei[d^<Pp] I -0y!Poi^)e'"^^^^ (A6) 

where 
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i 

For the case of a perfect lattice pQ^q) is 

Po{q) = Po{2nYY.5{q - K) (A8) 

Using ( |A8|) in (^) one gets formula (p.8|) . 

p(x)=podet[9„0^]Ee*^<^(^) (A9) 

Assuming that we are in the elastic limit daUp <^ 1 one can expand ( |A9| ) to get 

p{x) ^ po[l - d^uMi^)) + e^^(^-"(<^(^)'^))] (AlO) 



In (|A10|) one can replace u{(j){r, z), z) by ^(r, 2;) up to terms of order daUjS ^ 1. Note that 



in doing so u has negligible (suppressed by powers of a/^) Fourier components outside the 
Brillouin zone, and thus there is a complete decoupling between the gradient term and higher 
K terms. 

The same procedure can be carried on in the case where the equilibrium lattice of the Ri 
contains topological defects such as dislocations, vacancies etc. Suppose for instance that 
one wants to study a lattice with a fixed number of dislocations at prescribed positions in 
the internal coordinate of the lattice, i.e a network with a fixed topology (connectivity), but 
which is now allowed to fluctuate in the embedding space due to coupling to disorder and 
thermal noise. This is relevant for the physical situation of a flux lattice with quenched-in 
dislocations whenever one can neglect dislocation motion, i.e glide and climb. The equilib- 
rium positions Ri now correspond to a minimum of the elastic energy with the constraint of 
prescribed connectivity. The problem at hand is to analyze the extra small elastic displace- 
ments around the equilibrium position due to disorder. The density is still given by (|A6|) 
where po{q) is now the Fourier transform of the lattice of the Ri ([ATI). 

Upon coupling to disorder, the equivalent of the last term in (|2.13|) will be generated. 
After averaging over the random potential it reads: 

'^''^ ^'■cos{q{u''{x)-u\x)))S{q) (All) 



(27r)'^2T 

where S{q) = Po{q)po{—q) (with no averages) is the structure factor of the lattice without 
disorder, with fixed connectivity and in its equilibrium position. If in addition one wants to 
allow for topological defects to equilibrate one needs to perform some further average over 
connectivities, which is a difficult task. 



APPENDIX B: LINK WITH QUANTUM MODELS IN 1 + 1 DIMENSIONS 

The variational method can also be applied to study disordered one dimensional (space) 
interacting bosons or fermions. To see that fact we use a representation of operators in 
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terms of phase fields introduced by Haldane0. This representation maps the system into an 
elastic hamiltonian similar to the one used to describe flux lattices in ||. 
The single-particle creation operator is written: 

= [p(x)]^e'''(") (Bl) 

where p{x) is the particle density operator and 9{x) the phase of the \l/ field. To take into 
account the discrete nature of the particle density, one introduces an operator $ which 
increases by vr at each particle's location. The density operator then is 

1 d^(x) 

p(x) = exp[2im$(x)] (B2) 

^ m=-oo 

$ can be expressed in term of another operator </> by $ = irpox + (j), and po is the average 
density. The (j) and 6 fields obey the canonical commutation relations 

[(f){x),-e{x')]=i6{x-x') (B3) 
vr 



Using (|B1|) , the single particle operator becomes^ 

■q/\x) = [po + -V0(x)]i/2 e''"("''°"+^("))e*^(") (B4) 



For fermions the sum over m in ( p^ ) is only over odd m, whereas for boson the sum is only 
over even values of m. For fermions vrpo can be replaced in ( P4|) by kp, where kp is the 
Fermi momentum. The long- wavelength - low-energy properties of the interacting boson or 
fermion gas are described by the Hamiltonian@: 

H = ^l dx{{^){d,<pr + {vK){dJY} (B5) 

When going to the Lagrangian, one gets 

C=-^ I dxdr[\dr(t)f + (B6) 



which is obviously the same than the classical Hamiltonian ( p.2|) . 

From Galilean invariance one has {vK)/n = po/m, and v"^ = l/(fi:pom), where k is the 
compressibility, v and K are therefore functions of the interactions and incorporate all the 
interaction effects. The excited states of H are sound waves with phase velocity v. For the 
fermion problem the noninteracting case corresponds to v = Vp and K = 1. For repulsive 
interactions K < 1, whereas > 1 for attractive ones. For the boson problem — > oo 
when the repulsion between bosons goes to zero, and K decreases for increasing repulsive 
interactions. For the case of a 5 function repulsion, K varies from cx) to 1. K = 1 would 
correspond to an infinite on site repulsion. K < 1 can be obtained only if longer range 
interactions are considered. The coefficient K determines the asymptotic behavior of the 
correlation functionsS. For the bosons one gets 
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(^;,(r)^B(0)) = i?po(por)-i/(2i.^) 

(p(r)p(O)) = 2K,{2npor)-^ + ApHpor)-^''^ cos(27rpor) 



with some numerical constants A and 5. whereas the fermion correlation functions are 

(B8) 



(p(r)p(O)) = 2Kf{2TTpor)~^ + ApHpory^^'f cos(2A;^r) 



up to an angular part. The coupling to a random potential uncorrelated in space and time 
will give again a Lagrangian identical to ( 2.1ij| ). Terms of the form 



-^dis = E E / dxdrDp cos(2p(0'^(a;, r) - r))) (B9) 

where p is an integer and Dp some constants proportional to the disorder will appear. They 
correspond to Fourier components of the random potential close to 27rppQ. The Fourier 
components close to g = correspond to forward scattering and give terms similar to the 
gradient terms in ( p.l3| ). 

As discussed in Section ^ space and time uncorrelated disorder is relevant when K < 1. 
The transition point K = 1 corresponds to non-interacting fermions or equivalently to bosons 
interacting with infinite S repulsive potential. Disorder is thus always relevant for fermions 
with repulsive interactions, and always irrelevant for fermions with attractive interactions. 
Disorder is irrelevant for bosons with only finite (ie/ta-function repulsion, while it becomes 
relevant for bosons with sufficiently strong longer range interactions. One experimental 
realization corresponds to flux lines in superconductors, confined to a plane, which can be 
realized by proper alignement of the magnetic field. It was argued in Reli that this always 
lead to a glass phase, i.e K < 1. This does not seem correct to us. If the vortex line 
interaction was the sum of an infinite repulsive delta (i.e forbidden crossings) and an extra 
repulsive interaction of finite range A, then disorder would indeed always be relevant. But 
this is not the case, and the on-site interaction is finite, which presumably leaves room for 
both a glass phase {K < l)and a high temperature phase {K > 1) in this experimental 
system. 



APPENDIX C: ONE-STEP SOLUTION FOR D <2 



In this section we examine a one step replica symmetry broken solution. For simplicity 



we look at the model with a single harmonic (|3.11| ) We now search a one step symmetry 
breaking solution of the form 



a{v) = (To 
a{v) = ai 



[cr](f) =0 V < Vc 
[a]{v) = v^{ai - o-q) = 



(CI) 



Si 



V > Vr. 



With the form (pif) one gets 
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G{q) 
G{q,v) 



1 



Sia(g) 



v,{G^^{q) + j:i 
G{q,v) = croGl{q) v < v, 
[G]=0 V <Vc 



+ croGl{q) V > Vc 



{G-'{q) + Si 

To determine Vc one has to minimize the free energy 

T f d'^q 



(C2) 

(C3) 

(C4) 
(C5) 

(C6) 



F/{nmVt) 



2 J {2t:Y 
G(g,0) , 



cq'G{q) - \og{GM)) 



Gc{q) Jo V 
A /-i 



1 dv, 
— log 



Gcjq) - [G] 
Gc{q) 



(C7) 



where 



d'^q 



2mT Jo 



1 



B{0,v) = '^T J 
B{0,v) = 5(0, w > Vc) + 2T J 



{2'nYG-\q) + T., 
d'^q 



V > Vr 



Si 



{27rYG-^{q){G~^{q) + j:,)v, 



V < Vr 



To get the saddle point equations one has to differentiate (|C7|) with respect to o"o(g), Gc, 
Si and Vc- Differentiation with respect of CTo{q) gives G^^{q) = cq^. Differentiating with 
respect of 1/vc and Si and then replacing G~^{q) by its value one gets 



Si 



mT 



d'^q 
{2nY 



q 1 



log(l + ^) - 



(C9) 



Si 



cg^ cg^ + Si 



Note that ctq does not appear in these equations. One could in principle determine the value 
of (Jo by differentiating ( |C7|) with respect to G~^{q). However it is much easier to use the 
saddle point equations (|3.9| ). For d <2 one has 



^ a J (Ott 



9 L. 



mT 



ao = o-ie 



J (2ir)' 



(CIO) 
(Cll) 



Thus (Jo = for (i < 2. Note that in order to determine Vc one should not substitute in ( |C7| ) 
the expressions ( |C10|) before differentiating with respect to Si and 1/fc- Such expressions 
are only valid at the saddle point. 
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In (i = 2 the above equations read: 



Si 



A 



Si, 



TiTo'— log[l + -^ 



A 



(C12) 



where we used a circular cutoff A ~ c{27i/a)^. Let us denote by Tc = Attc/Kq and y = Si/A. 



(C13) 



We then have: 

T log(l+y) 
Tc y 

Al-T/T.) ^ ^gA l0g[l+l/] 

^ mrcAy(l + ?/)^/^= 

We are interested in the case ^ < 1 since ^ ^ a and thus there is clearly a one step 
solution for T < Tc. Note that for arbitrary A Si vanishes when T T~ and simultaneously 
Vc goes to 1. 

For d < 2 the above equations lead to 



d 



Si 



mT 



where 



Jd 



d'^q 



TT 



l-d/2 



{2tiY g2 + 1 2"^ sin(7rrf/2)r[rf/2] 



(C14) 
(C15) 

(C16) 



We can compare the free energy (|C7|) of the one step solution with the free energy of the 
symmetric one, that one obtains by taking f c = 1 and Si = 0. The free energy difference 



onestep 



- is 



^ = -(--1) 

mnil 2 Vc 



d'^q 
(27r)" 



Si 



cg2 + Si 



log 



.cg^ + Si 

9 



+ (1 



2mT 



-KlB(Q,Vr)/2 



(C17) 



where 5(0, Vc) is obtained from (|C8| ). Using the saddle point equations (p9|) , ( |C17|) simphfies 
into 



AF 

mnVt 



Si (1 



2Kl 



(C18) 



Thus whenever there exist a one step solution it has higher energy than the replica symmetric 
solution. 

Examining the equations |C14| one sees that indeed a one step solution survives in d = 1 
even while the replica symmetric solution is stable there. Note however that the one step 
solution changes nature and, as an approximation to the original problem, acquire some 
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unphysical features. In d = 1 it does not exist for small disorder at fixed T. When it 
exist the transition is discontinuous and Si jumps from zero to some non zero value. The 
equations are in fact very similar to a variational approach to roughening in d = 1 where it 
is obvious that the cosine term is irrelevant and the variational approach cannot be trusted 
beyond weak disorder, where it gives sensible results. Here it is also clear that q = Kq 
disorder (i.e the cosine term) is irrelevant at large scale for (i < 2, in the sense that it will 
not pin the d = 1 manifold. It will however renormalize the long-wavelength disorder which 
simply adds to the thermal roughening u ~ L^/^. The physical meaning of the one step 
solution for d = 1 being unclear, although it could be related to transitions in the dynamics, 
we will not considering further. Note finally that in ci = the one step solution does not 
exist at all ( since the mass is zero). 



APPENDIX D: VARIATIONAL CALCULATION ON THE NON-LOCAL MODEL 



If one does not use the decomposition ( p.9|) of the density, the variational equation for the 
off diagonal self energy iflU instead of (|3.2|) (in the presence of the extra disorder (|3.71|) ) 

(ya^hiqi) = -y At (- ^L)e-^s5^ cos(qM (Dl) 



As for the local model (|2.13|) one has to look for a non replica symmetric solution. Let us 



use the Poisson formula valid for arbitrary lattice of Ri. 

vol{R,) ^ 6{x - i?.) = E e'""" (D2) 

where vol{Ri) is the volume of the unit cell and the K are the reciprocal lattice vectors. 
One can then replace the discrete sum in ( pi| ) by an integral and insert the integration over 
X in formula (pil). One gets 



-tax...) = ^E/''"Ve--i^j^^(^ - _!l_)e-^oos(,,r) (D3) 

From this exact expression, BMY kept only the components with K = Q. As pointed out by 
BMY,ii'ii if one takes v small enough, then B{r = l,v) is large enough so that, for all the 
r contributing to (p^), one can replace B{r,v) by B{0,v). The idea is that B{r,v) varies 
much more slowly that r^, which is true in the elastic limit. Indeed one can easily see that 
for small v and large r one has B{r,v) oc log{v^/v) + Clog(r). Thus B{r,v) ~ implies 
that r ~ log(f^/f )^/^, and therefore B{r,v) ~ B{0,v). If one does so and since B(0,v) ^ 1, 
one can restrict oneself to = or i^' = ±Kq in (^^). The off diagonal part of the self 
energy becomes 



2 „, , 



2 



a(g, v) ~ A5(0, v)-'-^/^ + cig^e" V^(o.-) + C2e'^^(°'") (D4) 



The two first terms in ( P4| ) come from the K = part in ( P3|) . The term in A is due to 
the extra unphysical disorder and, as long as A is finite dominates the long range behavior 
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(g and large -8(0, v)). The other terms correspond to the same separation of the Fourier 
components of the random potential than in the discussion leading to our model ( p.l3| ). The 
g ~ components gives the same contribution than the terms in ( p.l3| ) in VMaVwf,. Indeed 
this K = term comes from 

( J drdr'6{r - r' + u^ir) - u\r'))) (D5) 

This term measures the smooth change in local density due to the slowly varying displace- 
ment field, and keeping only this term amounts to neglect the discreteness of the vortices in 
the original Hamiltonian. The third term in (^^) comes from the q ~ Kq component of the 
disorder, i.e. the part which has the periodicity of the equilibrium lattice. For the physical 
disorder A = ( p4| ) is identical to the one obtained from the local model ( |2.13| ). 

In fact neglecting the x dependence in B{x,v) is equivalent to the replacement = 
X — u{x), i.e identifying — u{(f){y)) with u{x) — u{y), which led to the model ( p. 131 ). 

Such a replacement becomes exact in the elastic limit ((Vm)^) ^ 1 which has been checked 
self consistently on our solution. 



APPENDIX E: RANDOM MANIFOLD IN D = 2 

Let us consider the original model ( p.l3| ) in d = 2 keeping all the harmonics. We will treat 
only the simpler case of isotropic elastic matrix, but for arbitrary m and lattice symmetry. 
Similarly to the study in > 2 we use the rescaled quantities ( p.27|) which satisfy the rescaled 



equations ( p.STp . In ci = 2 one has ^ = a^c^m/ATi'^/^Ko and = TxT/a?c. The equations 
read: 

[s]{y) = h{z) 

y = -h\z)/h{z) (El) 



where h{z) is defined in ( |3.3lD . It appears clearly on ( [El[ ) that y has a minimum value which 

^^°"V (E2) 



IS 

2 



27r / 



Thus the function [cr](f) must vanish for v < vq = v^yo = T/Tc which is identical to the 
value of the breakpoint Vc in ( |6.2| ) of the one step solution for the single cosine model. 
The asymptotic behavior at large distance will thus be identical to the one of the single 
cosine model studied in section [111 (J 1[ The transition temperature also turns out to be 



exactly the same = Attc/Kq . However, for y > yo, [s]{y) will increase continuously until 
the breakpoint y = yc above which it becomes constant again. Using equation (|3.22D the 
breakpoint = h{yc) satisfies 
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where we used the circular cutoff A = 271 /a. Using the equation (|E1|) one sees that Zc is 
determined as the solution of the equation Zc = f{zc) where we have defined the function 
f{z) as: 

Let us study the function f{z). Using the two hmiting forms for h{z), i.e 

h{z) = Ae""" z » 1 (E5) 

from p.39| and [3^^ , one finds that the function f{z) has three different behaviours depending 
on the values of z 

f{z)c^—^{az + \og{{i/af/A)) z»l 



fiz) ^ _(^_log(z) + hgii^aY/B)) {aUr^^^^^ « z « 1 (E7) 

■ttT 

f{z) ^ -^(e/a)2z(™+4)/2 ^ ^ (a/0^/("^+^) 

taking into account that ^/a » 1. 

The first regime corresponds to high temperatures. One finds no solution when T > 
and when T < Tc one finds Zc = nT \og{{^ / a)"^ / A) / ca^ {1 — T/Tc). In that regime i/c and yo 
are very close to each other and the solution is very similar to the one step solution of the 
single cosine model. There is no true random manifold regime. 

Solving the equation in the second regime one finds approximately 

IT I Tfi A~ 4 

ca^ 21og((4/a)7E) 

where we have used that in that regime y ~ (m + 4) / {2z) . The condition Zc 1 shows that 
this regime exists only at low temperature: 

T<T*= , , (E9) 

There all the harmonics contribute, j/c and i/q are very different from each other and there 
is a large random manifold regime, between the Larkin regime and the asymptotic regime. 
In this regime one has B{x) ~ A{x/^)'^'^ with u = v^m = 2/(4 + m). 

It is easy to see that the last regime of behaviour of f{z) never contributes. This is 
because one is restricted to values of z such that Vc = VcV^ < 1 which corresponds to 
Zc > ca^/(T7r). This regime where there is no breakpoint (/ ~ a), and therefore no Larkin 
regime, as already discussed after formula ( |3.56|) , arises at very low temperatures T < T'* ~ 
(caV7r)(a/0^/(^+™Vlog((e/a)V5) 
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